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I. INTRODUCTION 



A. BACKGROUND 



The nature or battle is rapidly changing. 



The development of modern warfare is outstriping our 
ability to manage it effectively. The speed, power and 
sophistication of modern weapons and the political and 
geographical complexities of potential battlefields in 
places like Europe and the Middle East have placed the 
commander in a rapidly moving, data-rich environment. 
[Ref. 1: p. 50] 



Many publication on the state of the modern battlefield 
use the same descriptive adjectives: intense, dynamic, 

mobile, rapidly changing, data-rich, etc. The quoted 

article goes on to stress that modern decision aids are 
needed to assist the battlefield decision maker . . by 

expanding his ability to rapidly and effectively analyze the 
data becoming available." [Ref. 1 : p. 50] However, to 

assist the force commander on the battlefield, much military 
research has emphasized the need to help the field commander 
and staff in deciding a course of action in critical battle- 
field situations [Ref. 2]. 

B. THE PROBLEM 



The commander’s decisions on the battlefield are prima- 
rily influenced by three factors: the environment, the enemy 

force, and the friendly force. The only factor that the 
force commander can directly control is the friendly force. 
If this control is to be effective, the commander must have 
intelligence about the environment, and the enemy force. In 
particular, the commander in the Airland battlefield needs a 
rapid, accurate assessment of the current enemy situation 
and the enemy's intentions. 
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The commander’s information requirements place an ever 
increasing responsibility on the intelligence analyst. Raw 
data from battlefield sensors must be converted to intelli- 
gence as rapidly as possible. Without a proper under- 
standing of the enemy's capabilities and intentions, the 
likelihood of selecting an appropriate course of action is 
extremely small. 

C. THESIS OBJECTIVES 

The primary objective of this thesis is to demonstrate 
the use of the Dempster-Shaf er theory of evidence on the 
battlefield to help correlate data for predicting enemy 
intentions. Specifically, the thesis will: 

1. Demonstrate the use of the Dem pster-Sha f er theory of 
evidence in a tactical military intelligence decision 
aid. 

2. Create a decision aid prototype that uses the 
Demps ter-Shaf er theory. 

3. Analyze techniques for reducing the computational 

complexity and calculation time required for 

Dempster-Shaf er. 

To understand the role of evidential reasoning on the 
battlefield. Chapter II of the thesis will discuss the 
tactical command and control process. It will also describe 
the problem of the intelligence analyst, unaided by an 

evidential decision aid, as he deals with the bombardment of 
battlefield data. Enhancement of the command and control 
cycle will be discussed, leading to an example showing the 
benefits of an improved cycle. 

In Chapter III, the concept or evidential reasoning will 
be introduced. Two methods £or„ evidential reasoning. 
Boolean Logic and Bayesian Inference, will be investigated 
to expose their shortcomings for use in battlefield 



13 



its 



situations. The Dempster- Shafer theory of evidence, 
advantages and disadvantages will be presented. 

Chapter IV will describe the Demps ter-Shaf er theory 
decision aid created to assist the military intelligence 
analyst. The aid functions in a cyclical process, receiving 
each bit of data from the user, combining it with previous 
data, and presenting its conclusions to the user. Reduction 
in computational complexity and calculation time of the 
Dempster-Skaf er method are also discussed. These computa- 
tional tecnnigues speed up the hierarchical summations 
required oy the Dempster-Shaf er theory, but they also reduce 
the scope of the problem. The reduction methods can be 
applied to evidential reasoning on the battlefield. 

Chapter V discusses the battlefield intelligence 
analyst's jab of Situation Development Analysis, a specific 
application area for the Dempster-S haf er decision aid. 'The 
enhancement of Situation Development Analysis by a 
Dempster-Shaf er decision aid will lead to benefits of 
improved command and control. 

A listing of the PASCAL code written for the decision 
aid in Chapter IV is contained in Appendix A. 
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II. BAT TLEFIELD CCHMAMD AND CONTROL 



The opportunity for battlefield applications of the 
Dempster-Shaf er theory of evidence stems from the need to 
shorten a commander's command and control cycle. With a 
shorter cycle, a force commander can recognize enemy inten- 
tions and strike quickly enough to disrupt the enemy opera- 
tion. With a quick reaction, the battlefield initiative can 

w 

be seized before the opposing force can react to the new 
situation. The now weakened and off balanced enemy force 
will be easier to defeat. This process of converting data 
into action more quicjtly than the opponent is referred to as 
turning within the enemy's decision cycle. 

Reducing the commander's command and control cycle time 
requires tue capability to react more guickly than ever 
before. The faster and deeper that the friendly force can 
interdict the enemy's forces, the more successful the 
o per at i on. 

Enhanced command and control has emerged as a solution 
for dealing with a dynamic and data enriched battlefield 
while reducing the commander's reaction time. This chapter 
will discuss the command and control cycle, the data flow in 
the cycle, current enhancements to the cycle, and present an 
example that shows the benefit of reducing cycle time by 
enhancing command and control. 

A. THE COaSAND AND CONTROL NETWORK 

Command and control is an extremely complex battlerieid 
function by which the commander and staff allocate 
resources, direct unit movement, and coordinate operations. 
Command, and control can be modelled as a continuing cycle 
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within a network of nodes as illustrated in Figure 2.1. 
Within this cycle are processes, depicted as the nodes, that 
can be enhanced to reduce cycle time. These nodes have been 
described in many ways by many experts £Ref. 3], but will be 
designated here as Collection, Interpretation, Decision, and 
Action. 
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Figure 2. 1 Command and Control Cycle Network. 

To understand the command and control network, and the 
cycling within it during battle, it is necessary to investi- 
gate each node of the network. Within each node, functions 
that can be enhanced or eliminated may be discovered. The 
improvement of node efficiency will reduce cycle time and 
help in the effort to interdict the enemy force as early as 
possible. 

The nodes of the cycle network depicted in Figure 2.1 
are as follows. 
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1. Collection 



The Collection node of this network describes the 
activity of gathering data about friendly and enemy forces 
by all means available. Miller and Cushman [Refs. 5,1], 
describe this activity as sensing, but sensing may include 
the recognition or interpretation of the data received. For 
purposes of this thesis, sensors only gather the data and do 
not modify it in any way. 

Collection is a continuing process and does not 
realize divisions between battles or battlefields. A sensor 
may receive data that may has no effect in its area of oper- 
ation. Yet, this same data may be of critical importance to 
the commander. Nonetheless, all of this data will be trans- 
mitted to the Interpretation node. 

2 • I nte rpre tation 

Interpretation is the combination, evaluation, and 
translation of raw data into intelligence. The result of 
this activity is an understanding of the battlefield situ- 
ation including enemy intentions. 

Interpretation, like Collection, is a continuing 
process. The intelligence analyst must deal with the 
continuous stream of relevant and irrelevant data arriving 
from the Collection node. Decision aids are used to help 

the intelligence analyst comprehend the entire battlefield 
picture presented by this onslaught of data. 

The arc connecting the Collection and Interpretation 
node indicates flow in both directions. As tne battlerield 
situation develops, sensors may be directed to change posi- 
tions or sense other areas of the battlefield. 

Consequently, the connecting are depicts coordination 
between these nodes. 
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3 . Decision 

After recognizing the enemy in 
force must react. This node describes 
sion process that chooses the course 
mizes possibility of friendly success 
course of enemy action. This success 
minimization of friendly losses, the 
losses, or by any measure of effective 
performance considered essential by th 

It is possible that the com 
enough information to justify a cours 
then reguest more data collection or 
other areas. The commander will inf 
node of his Essential Elements of 
Priority Intelligence Requirements (PI 
items of the utmost importance to the 
arc connecting the Interpretation a 
bi-dire ctional. 

4. Action 

The last node of the network r 
of troops into battle executing the 
the Decision node. This activity ma 
cated maneuver to a reactionary t 
encoun t er . 

The arc from the Decision to 
directional. This does not mean that 
dc not coordinate or respond to super 
the final orders ox the higher-level 
planning is accomplished. 

The command and control networ 
data from tne Decision node to the Ac 
2.1, and therefore, can only functi 
nodes have completed their processes. 



tentions, the friendly 
the commander's deci- 
of action which opti- 
given the interpreted 
can be measured by the 
maximization of enemy 
ness or any measure of 
e commander, 
mander does not have 
e of action. He may 
direct surveillance in 
orm the Interpretation 
Information (EEI) or 
R) , those intelligence 
force. Therefore, the 
nd Decision nodes is 



epresents the movement 
course of action from 
y range from a compli- 
actic in a surprise 

Action node is uni- 
subordinate commanders 
iors, but it depicts 
commander after all 

k ends in the flow of 
tion node, see Figure 
on after the preceding 
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Of course, the enemy must also cycle through a 
similar network during the battle. During operations, each 
force must constantly collect, interpret, decide, and act at 
all levels of command within the battlefield. The force 
that executes their cycle the fastest will have the battle- 
field advantage. 

E. ENHANCING THE NETHORK 

There are many areas within the command and control 
network that can be enhanced by some type of automation. 
Stewart, Boss, and Tiede [Ref. 4] have identified functional 
areas in any organization where the human processor can be 
improved by introducing automation. Many of these areas 

have applications in the command and control network. Tarle 

1 depicts these processes, the advantage gained by automa- 
tion, and the areas of the command and control Network where 
they apply. 

1 . Enhanc emen t Of Data Flow Through The Network 

The arcs which connect tne nodes of the command and 
ccntrol cycle network represent data flow. The Collection 
node passes data to the Interpretation node which passes 
data, or l n this frame, intelligence to the Decision node. 
Here, possible courses of action are analyzed to select the 
one that best fulfills the commanaer’s effectiveness 
measure. Next, the Decision node sends orders, which 

consist of the directions to units in the force, tc the 

Action node. These force units will then execute the chosen 
course of action. The force units then act or interact with 
the enemy forces creating more data to be gatnerea by the 
Collection node and the cycle continues. 

Much research has gone into the area of improving 
data handling and flow within the network [Ref. 5]. The 
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In the Interpretati 
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The Decision node is an area which has received 
substantial attention. In this node the commander and staff 
plan tuture courses of action, wargame these options, and 
then select the course with greatest probability of success. 
Enemy Courses of Action (ENCOA), Forces Comparison Model 
(FORCECOM), Contingency Screening Model (CONSCREEN) , and 
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Quick Screening Model (QUI CKSCREEN) are current models 
developed for the commander and staff to aid in this deci- 
sion process. With the abundance of models available in the 
military community related to decision making, this node was’ 
a good candidate for early enhancement by automation 
techniques. 



d. Action Enhancements 

Improving data flow in the command and control 
system, along with improving battlefield transportation will 
accelerate the purpose of the Action node ny getting troops 
to the required position on the battlefield in the most 
efficient manner as possible. 



C. THE CRITICAL NODE 

Which node in the command and control network holds the 
key to battlefield success? Clearly, the Collection node 
will have a very rich environment from which to gather data 
and transmit to the Interpretation node. Every unit on the 
battlefield will be forwarding data to be interpreted. 

The commander and staff may plan continuously for every 
conceivable enemy course of action. Forces will he posi- 
tioned to provide as much deterrent as possible to the 
e ne m y . 

Sub-units will plan for movement to each battle position 
for each friendly course of action directed by higher head- 
quarters. Therefore, recognition or interpretation of the 
enemy course of action becomes the limiting factor in the 
command and control cycle. 

If the initial enemy attack into NATO positions in 
Western Europe, or any activity” during war, can be 
predicted, the commander can commit his forces, especially 
the reserve forces kept for this purpose, to repel the 
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opponent. The Interpretation node is the trigger to this 
action. It must be accurate, effective, and as streamlined 
as possible. The entire network process hinges on the crit- 
ical role of the Interpretation node. 

How much of a difference will enhancement of this node 
by evidential reasoning, or any method for that matter, and 
the resultant reduction in command and control cycle time, 
make on the battlefield? 

D. BENEFIT OF SEDUCED COMA AND AND CONTROL CYCLE TINE , AN 

EXAMPLE 

An example [Ref, 5], is proposed to demonstrate the 
benefit achieved at a critical battlefield confrontation by 
reducing the command and control cycle time. 

In a European type scenario, units are deployed with the 
mission of maintaining some force ratio threshold (RTH ) in 
all sectors of operation in order to prevent an enemy break- 
through. Assume there is a predetermined critical sector in 
which the enemy has chosen to attempt such a breakthrough. 
The force Ratio in the Critical Sector (RC3) will be 
compute d as : 

RCS = R(t) / B (t) (2. 1) 

where ’R(t) is the total number of major t ank/anti-tank 
systems in the attacking (Red) force at time (t), and B (t) 
is the same identifier for the defending (Blue) force. Both 
force describers are functions of time to allow for attri- 
tion and reinforcement by reserve forces. The Red forces 
are divided between critical and non-critical sectors. The 
Blue forces have in position a reserve force containing part 
(percentage of force = x) of the total force. See Figure 
2.2 for force dispositions. 
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Critical Sector 



R(t) 

(1 -x) «B ( t) 

Reserve Force 
x* ( E (t ) ) 

* 

R ( t) = Attacker Strength at Time (t) 
B ( t) = Defender Strength at Time (t) 
x = Percentage of Force in Pveserve 



Figure 2.2 Battlefield Layout. 

If RTH is the threshold ratio necessary to prevent .a 
breakthrough in the critical sector, the defender must keep 
RCS < RTH. To accomplish this standoff, the command and 
control cycle must function in a timely manner so that 
reserve forces are committed at the correct moment in time 
and the required threshold ratio is maintained. 

Figure 2.3 shows a hypothetical time comparison of the 
attacker/defender schedules for respective command and 
control cycles. Figure 2.4 shows hypothetical linear force 
ratios indicating mission deficiency on the part of the 
defending force due to the surprise gained by the attacker’s 
commitment in the critical sector. 

Reduction of the defender’s command and control cycle 
will increase B (t) , increase attacker attrition, and reduce 
the mission deficiency, see Figure 2.5. The time sequence 
used in Figures 2.4 and Figure 2.5 correspond to those 
schedules depicted m Figure 2.3. 



25 



Attacker Activities: 




I 



I 



Tl-Attacker decides 

T2-Movement of main 

T3-Lead elements of 
critical sector 

T4-Last elements of 
critical secto r 



to commit main attack force 

attack force begins 

main attack force arrive in 

main attack force close in 



Attacke r -- | | | | 

Activity 

T1 T2 T3 T4 



Collect Interpret Decide Act 

Defender | | | | | | 

Activity 

tl t2 t3 t4 t5 t6 



Defender Activities: 

tl-Defender detects movement of main attack force 

t2-Defender interprets time and location of main 
attack 

t3-Defender decides to commit reserve force to 
critical secto r 

t4-Lead elements of reserve force arrive in 
critical sector 

t5-Last elements of reserve force . close in the 
c r i t i ca I sector 

t6-Reserve fully deployed into new defensive 
positions in the critical sector 



Figure 2.3 Time Comparison of Attacker 
and Defender Cycles ( notional) . 



Although this is a very simplistic and hypothetical 
situation, it nonetheless demonstrates the advantages that 
occur when the command and control cycle is shortened Dy any 
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node in the network. The earlier the critical sector can be 
reinforced, or the deeper the attacker can be interdicted, 
the more likely the ratio in the sector can be reduced to 
RTH and a breakthrough prevented. Consequently, any reduc- 
tion in cycle time can be directly equated to an increase in 
effectiveness of the force. Therefore, reduction of the 
Interpretation node's processing time by a decision aid 
using evidence combination will have a positive effect on 
the battle. 

To further reduce cycle time, it is possible to 
initially have a contingency plan allowing for a major 
attack in all sectors of the front. Once the critical 
sector is identified by the Interpret node, the plan for 
that sector would be executed. This contingency would 
almost remove the Decision node from the command and control 
network. This planning concept would reduce the cycle time 
by that former amount allocated to the Decision node and 
further reduce the mission deficiency. 
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Figure 2.4 Comparison of Attacker and Defender Force Ratios. 
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III. EVIDENTIAL REASONING 



A. BACKGROUNE 

As battlefield information flows into the Interpretation 
node of the Command and Control cycle, it must be combined 
or fused with prior information to update the understanding 
of the enemy situation. As previously discussed, this 
process is currently being accomplished by a man-in-the-loop 
system that can slow down the entire command and control 
cycle [Ref. 1]. . Methods of evidence combination that can 
enhance, not replace, the human inference process will now 
be discussed. 

The nattlefield situation of Chapter II described tne 
commander's desire to recognize the sector of main attack. 
The determination of the enemy's intentions can be viewed as 
a test of hypotheses consisting of all mutually exclusive 
sector combinations. Given the knowledge that an attack is 
imminent, such as the first battle in a European scenario, 
it is the analyst's job to accept or reject these hypoth- 
eses. To accomplish this task, the analyst must place 

values on each likely sector, or hypothesis, indicating the 
probability of an attack in that sector. 

The comparison here of battlefield reasoning to 

hypothesis testing is logical. The analyst has a set of 
hypotheses, composed of tne sectors of possible attacks, and 
their multiple conjunctions. These hypotheses indicate 

attacks over any one, or any combination, of the sectors in 
the force commander's zone of responsirili ty. 

An example of this bat tlefield- s ituation will now be 
described. It will demonstrate the use of an evidential 
reasoning process in the Interpretation node. This example 
will be used throughout the chapter. 
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To begin, suppose the unit using this reasoning enhance- 
ment is a U.S. Division. The division's area of responsi- 
bility is divided into three brigade sectors (1,2,3). The 
enemy attack, assuming this is the state of the battle, 

could occur in Sector 1, 2, or 3. If it occurs near a 
brigade boundary, or is wide enough to cover more than one 
sector, the attack could occur in combinations (1,2) or 
(2,3). If a divided attack occurs, or attacks from adjacent 
enemy forces occur, then (1,3) is a possibility. 

The battlefield, or population, is sampled via the 
unit's sensors, and the sampled evidence leads the analyst 
to accept or reject the hypotheses of attack locations. In 
the case of battlefield sampling, the sampling process takes 
place over time. 

This chapter will discuss the quality of the evidence 
presented to the decision maker, the use of Boolean and 
Bayesian methods to evaluate the analyst's hypotheses, and 
an in depth look at a technique for evidential reasoning 
based on the De mpster-Shafer theory of evidence. 

B. THE QUALITY OF EVIDENCE 

In the tactical environment, the Interpretation node 
will receive information from many sensors, both human and 
machine. This information is inherently uncertain, incom- 
plete, and sometimes inaccurate. Although less than 

optimal, this situation is the nature of the battlefield and 
the nature of the evidence received. The terra evidence 
becomes appropriate here since the information received will 
be a basis fcr conclusions or judgements, not a clear answer 
to any one hypothesis. 
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1 . Unc er tain E vide nce 

Battlefield sensors cannot describe their sample in 
precise detail and, therefore, create inherent uncertainty 
in the reports they generate. Some sensors can indicate 
movement direction, size of the element, of type of unit 
observed. Other sensors may only report that a unit is 
moving and may not detect the size or exact location of the 
target. None can give a complete description of the event. 
The sensor or operator is always uncertain of many attri- 
butes of the target. As lowrance and Garvey recognize 
[ Bef . 6], the evidence tends to lend varying degrees of 
support to one or more hypotheses rather than completely 
specify the event. 

2. Incomplete Ev i de nce 

The sensor information will also be incomplete. The 
battlefield sensor can only "view" its assigned sector of 
search. It can describe what it sees, but cannot lend 
evidence - to what it cannot see on other parts of the battle- 
field. The analyst should realize that this incompleteness 

exists, and direct movement of sensors or change sensor 
search areas to receive a more complete battlefield picture. 

3. Incorrect Evidence 



The third characteristic of the information 
collected by the sensor is that it could be incorrect. The 
operator, interpreter, or soldier reporting could be 
completely mistaken in their spotting, or the enemy could be 
using deceptive techniques to confuse the opposing force. 

For these three reasons: uncertainty, incomplete- 

ness, and i nccrrectness, the hypotheses or propositions can 
only be attributed degrees of support based on the evidence 
received. No one piece of information can be accepted as 
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complete truth. Therefore, it is the iterative process of 
combining information from all sources through time, that 
will lessen the damaging effects of poor-quality evidence 
and produce good intelligence. To accept or reject the 

propositions based on a single, cloudy piece of evidence 

would certainly bias the entire intelligence prediction 
effort . 

What statistical methods can be used to evaluate the 
evidence in terms of lending support to the acceptance or 
rejection of the hypotheses? Three methods of evaluation 
will be discussed. These are: Boolean Logic, Bayesian 

Inference, and Evidential Reasoning using the 

Dempster-Shaf er theory of evidence combination. 

C. TWO COSMOS HETHODS FOB EVIDENCE EVALUATION 

1 • Boolean Logic 

In Boolean logic the hypotheses or propositions can 
only be represented as True or False. Varying degrees of 
support are not accepted, and any information relative to 
the hypothesis would have to be interpreted as total support 
cr total negation. 

As Lowrance and Garvey indicate [Ref. 6], Eoolean 
logic cannot capture the partial belief in hypotheses gener- 
ated by the coarse evidence received. The battlefield will 
never be an area for clear cut decisions in black and white, 
but will always tend towards decisions that deal with the 
" grayness" of the evidence. 

Continuing with the Division example, a report is 
received of a small unit, an enemy motorized rifle company, 
moving towards the Forward Line of Troops (FLOT) of 
Sector 1. This activity could be the advance of the enemy's 
reserve force indicating a breakthrough attempt, or it cculd 
be a feigning action, cr only a partial repositioning of 
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troops. Based on this activity, the analyst cannot, 

100% certainty, predict that the main attack is coming in 
Sector 1 (True) or that it does not indicate an attack in 

Sector 1 (False) . Nor does the evidence negate or support 
the possibility of attack in the other sectors. 

2. Bayesian I nfere nce 

Bayesian methods have been proposed as a basis for 
several decision aids [Refs. 7,8]. However, these methods 
have inherent ineff iciencies in dealing with disbelief, 
information supporting the compliment of the hypothesis, and 
a priori probabilities. 

a. Bayesian Formulation 

A Bayesian approach [Ref. 7] would consider 
various hypotheses, such as an attack in Sector 1 (SI) given 
some datum (D) from a battlefield sensor. The probability 
that (SI) is true given the data would be: 

P (5 1 | D) = P (D | SI ) • ?(S1) / P (D) (3.1) ' 

where: 



P (SI J D) = posterior probability of the (3.2) 

hypothesis given the ooserved datum. 

P (D | S 1) = probability of the datum (3.3) 

given the hypothesis. 

P (51 ) •= prior probability of the hypothesis (3.4) 

before the datum. 

P (D) = probability of the datum occurring. (3.5) 
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The task of determining just this last prob- 
ability of a datum occuring on the battlefield, is more than 
formidable. Given a fixed set of target describers : 

TD = (tfi ( 1) ,td (2) ,td(n)) (3.6) 

Where td(1) and td (2) could be map grid coordinates in 
easting and northing, td (3) could be the type of unit, td(4) 
the size of the unit, td (5) the direction of movement, etc. 
The analyst would then express partial beliefs over TD by 
distributing belief to the elements of TD. 

for example, a report is received about a tank 
company in the intelligence analyst's area of responsi- 
bility. He must now determine the probability of detecting 
the tank company at the reported location. In this case TD 
= (td ( 1 ) , td (2 ) ) , the set of grid coordinates in the division 
zone. Probabilities must be mapped to all grid coordinates 
in his area which are maneuverable by tanks. The analyst 
then adds the probabilities for td(1) and td (2) corre- 
sponding to the reported grid coordinates of the tank 
company to determine the probability of the datum occurring. 

Probabilisticly, a piece of evidence will map 
the propositions in TD to the closed interval [0,1]: 

m : TD — > '[0,1] (3.7) 

w her e 



SUM m (td) = 1, for td in TD 



(3.8) 



In other words, the conjunction of all evidence from the 
mapping must equal one, the basis for a probability 
statement. 
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Then for any proposition defined over TD, such 
as the report of the tank company (TANK) , the probability of 
occurrence is: 

for all td in TD , (3.9) 

Prob ( TANK ) = SUM m(td), td in (TANK) 



The probability of proposition (TANK) is determined by the 
sum of the probability of all location possimiities that 
are elements of the proposition. 

It follows that: 



Prob (TANK) = 1 - Prob (not TANK) 
and since the environment sums to 1: 



(3. 10) 



SUM m (td) = 1 , td in TD 



(3. 11) 



Thus all probability not in (TANK) would lie elsewhere in 
TD, as seen in Equation 3. 10. The inherent problem with 
this approach is that the sensor operator, the intelligence 
analyst, or an expert, must determine each a priori prob- 
ability for the partitioning of TD into its elements 
{ td (.1 ) , td (2 ) , . . . td (n) ) given by the mapping m. 

This mapping would not be a great problem given 
a rich data base for a well-defined environment. However, 
on the battlefield, the sensors will be receiving data on 
micro events that may on.i.y occur once, and depend on time, 
weatner, terrain, or any other target descriptors used in 
(TD) . To compute a priori probabilities given the general 
nature of the situation would be an endless task and may not 
be acceptable for time-critical tactical decision making. 
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Using equation 3.1 and this chapter's example of 
a main attack against a Division, consider the Bayesian task 
of determining the probability of an occurrence.. If datum 
(D) is the report of the advancing motorized rifle company, 
then P ( D) is the probability of a such a company advancing. 
Immediately the analyst is in a predicament. He could make 
a calculated guess about this value if an extensive data 
base existed. However, there is no data base, and the 
analyst definitely does not have the time and may not have 
the expertise to concern himself with these detailed 
parameters. 

Further, what specific type of company was 
observed by the sensor? If it is a reconnaissance company, 
then the protability of this datum occurring could be high. 
Next, the probability of attack given this datum (P(S1|D)) 
must be considered. There are now more factors and prob- 
abilities with which the analyst must concern himself, most 
of which are not known. 

If the advancing company is a second echelon 
element, then there would be a low probability of the unit, 
being in this forward area, unless i t is the advance of a 
breakthrough attempt. This probabili Stic predicament could 
go on and on. Because of the many unknowns, Bayesian infer- 
ence may not be the most desirable method with whicn to deal 
with battlefield evidence combination. 

Also, a Bayesian supporter would say that if 
evidence supported two mutually exclusive propositions, and 
there was no reason to consider either over the other, each 
should be assigned equal probabilities. So, if evidence 
supports proposition (X or Y or Z) , with probability 0.6, it 
supports individual propositions (X) , (Y) , (Z) with prob- 
ability 0.2. As a result, there is a twofold support of the 
disjunction of any two of these propositions over the other. 
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lead to the theory of evidence proposed by Artnur Dempster 
and Glenn Shafer. This is not to say that Bayesian 
Inference should not be used on the battlefield. On the 
contrary, if the probabilities for the use of Bayesian 
Inference are available and if the analyst can distribute 
probabilities to single elements, then definitely Bayesian 
Inference should be used. However, if these conditions are 
not met, and tremendous assumptions would be required to 
meet the Bayesian prereguisi tes, then Dempster- Shafer should 
he considered as an alternative. 

D. THE DEHPSTEB-S HAFEB THEORY 

This theory of evidence comoination was conceived by 
Arthur Dempster [Ref. 9], and later developed by one of 
Dempster's students, Glenn Shafer [Ref. 10]. It is a theory 
of evidence because it deals with support of propositions 
based on evidence. It allows for quantifying ignorance, or 
lack of knowledge, as well as uncertainty. It uses a term, 
plausibility, to indicate lack of belief in a proposition 
rather than suggesting support of the compliment of the 
proposition. The term. Belief, is used to indicate support 
for any proposition. 

1 * For m ula tion 

The Frame of Discernment, the set of ail mutually 
exclusive propositions, is represented by: THETA. The 

domain of THETA is the set of ail possible subsets of THETA. 

Domain Size = 2 exp |THETA| subsets (3.12) 

(2 raised to the magnitude of THETA subsets) 

An example of a Frame or Discernment would be the 
three brigade sectors of the division zone in which the main 
enemy attack could occur. In this case THETA = (1,2,3), 

assuming 3 brigades in a division. 
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The domain would be : 



(1,2,3) (1,2) (1,3) (2,3) (1) , (2), (3) , () 

which are all the subsets of THETA. Although the Null Set, 
(), will always be a subset of THETA, it will never have 
Belief and will not be listed as a subset in the rest of the 
calculations in the chapter. 

The mapping of probability assignments to the propo- 
sitions is acne by a Basic Probability Assignment (BPA) , 
referred to as Mass (M), which satisfies: 

M (Null Set) =0 (3.13) 

and 



SUM M (i) = 1, i in THETA (3.14) 

(the sum of Masses over THETA = 1) . 



Support *for any proposition (X) is given by 
Eelief(X), or Bel(X), defined as: 

Bel (X ) = SUM 1 (i) , i in X (3.15) 

Eelief is the sum of all the Masses of all subsets of the 
proposition. The Belief of (X) is the measure of the Mass 

constrained to stay somewhere in (X) [Ref. 12]. Belief of 
an attack, in sectors (1 or 2 ) would egual: 

Bel (1,2) = M (1,2) + M(1) + M (2) 

Related to Belief is Plausibility or the Upper 
Probability Function defined as the total probability Mass 
that has potential for moving into (X). 



40 



Plaus (X) 



1 - Bel {not X) 



(3. 16) 



In the example, 

Plaus (1,2) = 1-Bel[not (1, 2) ] = 1-Bel(3) 

The Doubt of (X) or Dou(X) : 

Dou(X) = Bel (not X) (3. 17) 

is the measure of probability forced to stay out of A. 

Dou ( 1, 2) = Bel[ not (1,2) ] = Bel (3) 

It follows that; 

Plaus (A) > Eel (A) > M ( A ) (3.18) 

In contrast with Bayesian Inference, Dempster-Shaf er 
allows motion of Masses throughout the frame or discernment 
since each Mass need not be constrained to single elements 
within IH5TA. Therefore, no requirement exists to commit 
Masses to elements past tne level of recognition contained 
in the report, constrained by the sensor’s ii mite i 
capabilities. 

If evidence received indicates an attack in Sector 1 
or Sector 2, the support need not be divided between the two 
proposi tions , (1) and (2) . If movement of enemy forces 

towards the front occurs on a road bisecting the two zones, 
it is not necessary to say that; 

P (1) = P (2) = 1/2 • P (1,2) 

but instead, the evidence can be assigned to the superset: 

( 1 , 2 ). 

By using the two values of Belief and Plausibility, 
support for a proposition or hypothesis can me expressed by 
an interval as follows. 
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Evidential Interval (El) = [ Bel (X) , P la us (X) ] 



(3. 19) 



where the difference, Plaus(X) - Bel(X), can be referred to 
as the Ignorance remaining about (X) : 

Ig (X) = Bel (X) - Plaus(X) (3.20) 

If Ig (X) = 0, there exists no Mass available to move 
into (X). Further, if Ig(X) = 0 for all propositions, the 
system is Bayesian. This is true since this would require 
that all Masses be distributed to singletons in the frame of 
THETA. For this reason, Bayesian Inference can be described 
as a subclass of the theory cf Belief functions [Ref. 14]. 
Table 2 show some examples of these Evidential Interval 
v al u es . 



•TABLE 2 | 

EXAMPLES CF DEMPSTER- SHAFER EVIDENCE INTERVALS 



X (0,1) 


=> 


_ , 1 

No knowledge at all about (X). j 


X (0,0) 


= > 


(X) is false. 


X (1,1) 


=> 


(X) is true. 


X (0.25,1 ) 


= > 


Evidence provides partial 
support for (X) . 


X (0,0.85) 


= > 


Evidence provides partial 
support for (not X) . 


X (0 .25,0. 


85) => 


Probability of (X) is between 
0.25 and 0.85;i.e., tne evidence 
simultaneously provides support 
for both _(X ) and (not X). 
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In the continuing Division example, a report is 
received indicating a strong possibility of attack in Sector 
(1 or 2) represented as (1,2) and a slight possibility of 
attack in Sector 3, represented as (3) . The analyst assigns 
Mass (1,2) =0.6 and Mass (3) = 0.2. The remainder of Mass 

(0.2) cannot be assigned elsewhere, so is assigned to the 
set representing the entire frame, Mass(1,2,3) = 0.2. These 

Mass values must be determined in some way by the analyst. 
He may use probabilities derived from prior analysis of 
enemy tactics if available, or may use nis own judgement 
based on knowledge of enemy tactics. The latter method will 
be further discussed in Chapter V. For now assume that the 
analyst has determined these values. 

Belief is now limited to those sets whose subsets 
have Mass or: 

Bel (1,2, 3) = M ( 1 , 2, 3) + £1(1,2) + M(1,3) + M(2,3) 

+ M ( 1 ) + M ( 2) + M ( 3 ) 

= 0.2 + 0.6 + 0.0 + 0.0 + 0.0 + 0.0 + 0.2 = 1.0 

Bel (1,2) = M ( 1, 2) + M ( 1) + M(2) = 0.6+ 0.0 + 0.0 = 0.6 

Bel (2,3) = M (2,3) + M ( 2) + M(3) = 0 + 0 + 0.2 = 0.2 

Bel (1,3) = i 1 (1,3) + M(1) + M (3) = 0.0 + 0.0 + 0.2 = 0.2 

Bel (3) = M (3) = 0. 2 

Bei(1) = Bel (2) = 0.0 (They have no Mass assigned) 

As seen in the above example, since THETA exhausts all 
proposition possibilities, the Belief in THETA is always 
equal to 1.0. The Plausibility and Evidential Interval for 
(1,2) w ould be: 

Pla us (1,2) = 1 - Bel (not ( 1 , 2) ) = 1 - Bel (3) 

= 1 - 0.2 = 0.8 

Evidential Interval = [ Bel ( 1 , 2) , Pla us (1 , 2) ] 

= [ 0. 6 , 0 .8 J 
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The Ignorance would be: 



Ig (1 ,2) = 0. 8 - 0.6 = 0.2 

The Evidential Intervals (El) for all the sets are: 



El (1 ,2) = [ 0.6, 0.8 ] 

El (2,3) = [0.2, 1.0] 

El (1 ,3) = [0.2, 1. 0 ] 

El ( 1 ) = [ 0.0, 0.8] 

El (2 ) = [ 0 .0,0. 8 ] 



El (3) = [ 0.2, 0. 4 ] 



2. Combination of Evidence 

Dempster-Shaf er allows for the combination of 
evidence from knowledge sources or sensors on the battle- 
field. Given two Mass assignments Ml and M2 new Masses are 
computed by the orthogonal sum Ml a) M2 (where a) represents 
orthogonal sum) defined as: 

M (Null Set) = 0. 0 (3. 21) 

M ( A) = {SUM [ M 1 (X) • M2 (Y) ]} /K (3.22) 

where (X & Y) = (A) ; (S represents intersection) 

K = 1 - SUM [Ml (X) • M2(Y) ], (X S i') = Null Set 

= SUM [M1(X) • M2(Y) ], (X S Y) * Null Set 

From Eguation 3.21, the Mass assigned to the null set must 
equal 0. This is accomplished by the normalization of the 
Masses assigned to all other sets. The normalization 
factor, K, is equal to 1 - the Mass -assigned to the null set 
after combination. K is also equal to the sum of the Masses 
assigned to the' subsets of THETA (less the Null Set) after 
combination . 
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'The new Mass / or combined Mass assigned to (A) , is 
the orthogonal sum of the Masses (divided by the 
normalization factor) where (X & Y) = (A) . This orthogonal 

sum combination technique can be further shown by the unit 
square in Figure 3.1. 



M» (A) 




Mi (A) x M? (X ) 
dedicated to 

A tnt x 



l 

Mi (B) 



Mi (C) 




Unit Square 



Figure 3.1 Unit Sguare Example of Orthogonal Sum. 

Continuing with the example, a second report is 
received by the division intelligence analyst for which he 
assigns Masses as follows: 

M2(1,3) = 0.3 
M2 (2) = 0.3 
M2 (1 ,2,3) =0.4 

(where M2 designates second report). 

Figure 3.2 shews the orthogonal sum of these example Masses. 
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Figure 3.2 Orthogonal Sum of Example Masses. 



The Normalizing Factor (K) in this example is: 

K = (1-0. 06) = .94 

(1- Mass assigned to the Null Set) 



New normalized Masses are: 



M (1, 2, 3) = . 08/. 94 

M (1 , 2) = . 2 4/. 9 4 = 

M (1 , 3) = . 06/. 94 = 



. 085 



255 

064 
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3(2,3) = 0-0 (not an intersection set 
in orthogonal sum) 

M (1 ) = . 18/. 94 = . 192 

3(2) = (. 06 + . 18)/. 94 = .255 

3(3) = (. 08 + .06)/. 94 = . 149 

New Beliefs are computed as follows; 

Bel ( 1,2,3) = .085 + .255 + .064 + . 192 + .255 + .149 
= 1.0 

Bel ( 1, 2) = . 255 + . 192 + .255 =. 502 

Bel (1,3) = .064 + . 192 + . 149 = .40 5 

Bel (2,3) = 0.0 + .255 + . 149 = . 404 (has belief even 
though nc Bass, due to Mass in its subsets) 

Bel ( 1) = .192 
Bel (2) = . 25 5 
Bel (3) = .149 

The new Evidential Intervals are: 

El ( 1,2) = £.502, .851 ] 

EI(2,3) = £. 404 , .808] 

El ( 1 , 3) = £. 405, .745] 

El (1 ) = £ . 192 ,.596 ] 

El ( 2) = £ . 255, .595] 

El (3) = [ . 149, .498 ] 

These results show strongest Belief in the pair 
(1,2) due mainly to the assignment of Mass = . 6 to this set 
from Ml. Belief in (1,3) and (2,3) are approximately equal, 
.405 and .404 respectively. Belief for single elements is 
strongest for (2) due mainly to its Mass assignment from Ml. 
The final Beliefs follow from the Mass assignments, which is 
appealing. 

Also note the changes in the Evidential Intervals 
after combination. The intervals have narrowed for most 
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sets and no set has a Plausibility of 1-0. The Ignorance 
(Plausibility - Belief) for all sets except (1,2) has 
decreased. 

3. I nd ep en dence of Knowl edg e Sources 

A point of interest in Dempster's rule of combina- 
tion is the independence of knowledge sources or indepen- 
dence of reports from the same source. How does this method 
deal with multiple sightings from same or different sensors? 

If the sighting or sensor data come from different 
knowledge sources, then this evidence can be considered 
collected ircn independent elements- Sensors will occupy 
different terrain positions and have varying operating char- 
acteristics and capabilities. As such, these knowledge 
sources will derive enough independence for tnis method to 
produce desired results. 

If the same sensor reports the same data, these 
reports are net independent in the sense that they come from 
the same source. However, if the same sensor reports on the 
same unit of activity, but the location or any other repor- 
table characteristic of the unit changes, then this is 
considered sufficiently independent. An example would be 
reports on movement direction of a target. This reporting 
would be considered independent for the purpose of indi- 
cating a conndence towards or away form a hypothesis 
[Bef. 11 ]. 

If the same sensor is reporting on an activity that 
does not move, a higher headguarters location, then this 
should ne used as new evidence. It confirms headquarter 
location and also the fact they are not moving (shown by 
many reports from the same sensor) and therefore may indi- 
cate something to the analyst to be included in the evidence 
combina tion . 
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4. Some Points of Contention with De m pste r - Shafer 

Conflict over the Dempster- Shafer theory of evidence 
arises from the following: 

a) Normalization. This is the normalizing of resultant 
Masses after combination. Mass that would go to the 
null set is ignored and Masses going to THETA are 
normalized . 

b) Total Conflict. • When total conflict of evidence 
’.occurs. Mass for all sets sums to zero. The normal- 
izing factor (K) is egual to zero, and the attempt to 
divide by K, of course, fails. This case is the same 
as sending all Mass into the null set. 

However, there are several explanations to reduce 
concern over these events. 

(1) Normali za tion. Consider Masses tending to 
support conflicting propositions, propositions that Lave no 
intersection. The Mass in this event flows into the null 
set indicating conflict between Masses from Knowledge 
sources. Since this conflict occurs due to uncertainty 
about the situation at hand, as more evidence is received 
conflict will diminish. 

Also, as certainty toward the correct 
proposition increases so does ability to decrease the number 
of elements in THETA. Then, less conflict will occur in the 
evidence from sensors and the need for an evidence combina- 
tion technique that deals with uncertainty diminishes. But, 
this is not the case on the battlefield, where the sensors 
are spread over great distances across division fronts. 
These sensors are directed toward different areas of the 
front and a great diversion 'of information is desired by the 
intelligence analyst through a wide variety of contacts. 

Normalization, then, is just a means of 
dealing with the conflicting nature of evidence. The 



49 



measure or weight of conflict is represented by the magni- 
tude of the normalizing factor (K) . K indicates how much 
conflict occurs between the current Mass assignment and the 
resultant Mass values of previous assignments. If the set 
of all possibilities is the Frame of Discernment, the null 
set cannot occur. The Mass assigned to it represents 
conflict, a normal occurrence when so many sensors report so 
varied data. The Mass cannot remain in the null set, there- 
fore all Masses are normalized. 

A large measure of conflict occurs in the 

following example: 

Ml (A ) = .99 
M 1 ( 3 ) = .01 
M2(C) = .99 
M2(B) = .01 

M 1 s> M 2 yield : 

Ml (B) = .0001 
Ml (Null Set) = .9999 

Normalized results are: 

X = 1 - .9993 = .0001 

Ml o) M2 (B) = .0001/. 000 1 = 1.0, a questionable result 

However, if in fact (A,B,C) were the only 
possible results, this conclusion is logical. Only through 
comoining this oattiei ieid-q uali ty evidence and dealing with 
inherent uncertainty can the analyst reach conclusions about 
the hypotheses. If the rigorous methods lixe Bayesian 
Inference cannot be used, then Dempster-Shaf er seems to 
offer a logical alternative. 

(2) T otal Conflict. The second shortco ming, 
total conflict of evidence, is caused by total combined Mass 
going to the null set aborting any combination effort. This 
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occurs by assigning total Mass from two knowledge sources to 
sets with no intersection, such as: 

Ml (A ,B) = 1.0 
M2(C,D) = 1.0 

Ml 5) M2 = = > K = 0 

since : 

[ (A, B) & (C,D) ] = Null Set 

However, is it possible to be 1003 certain 
of (A, B), yet at tne same time, to be 100?f. certain of (C,D) ? 
No, this situation and resulting conflict are unacceptable 
and can only be resolved by the realization tnat, again, the 
analyst is dealing with uncertainty. Inaccuracies in knowl- 
edge sources cannot allow assignment of probability Masses 
in this manner. 

As further proof, consider the following 
example [Ref. 12]. A fair die is rolled and knowledge 
source one places all Mass in the proposition that the 
number is even. Ml (EVEN) = 1.0. Conflicting evidence from 
knowledge source two places all Mass in the proposition that 
the number is odd, M2 (ODD) = 1.0. The result of combination 
assigns all Mass to the null set since (EVEN & ODD) = Null 
Set, and the combination fails. 

This example shows that this occurrence 
would tend to violate the assumption of uncertainty of the 
evidence and also falsify any logical reasoning process. 
The analyst cannot let this occur under any method of 
evidence combination. 

E. BAYES* RUIE OF CONDITIONING 

In a battlefield situation, using Bayes* conditioning 
rule to combine evidence does not seem to generate more 
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satisfying results than those achieved by the 
Dempster-Shaf er method. 

A probability distribution m would be transformed by 
Eayes' Rule to m ' by the receipt of additional information. 
It is necessary to restrict the domain of m to elements of 
(X) when using this rule. Bayes' rule is: 

for all i in THETA: (3. 23) 

m'(i) = 0.0, if i not in (X) 

m ' (i) = m (i) / (1-k) , if i in (X) 

where: 

k = SUd m (i) < 1, i in (X) (3.24) 

Or from equation 3.1: 

P(S1|D) = P ( D | S 1 ) • P(S1)/P(D) 

Using the results of the prior knowledge source report 
outcome, P(S1), P(S1|D) is updated by the probability of the 
new datum, P(D), and probability of the datum given the 
hypothesis, P(D|S1). 

A more tractable form of the eguation can be constructed 
using likelihood ratios where: 

P (SI ID) /P (£2 J D) = [P(S1)/P(S2) ]• P (D | S 1) /? (D | S2) or 
P (SI |D)/P (S2 | D) = [P(S1)/P(S2) ]«L(D|S1:S2) 

where L(D|S1:S2) = P (D ] S 1 ) /P (D| S2) and is called the likeli- 
hood ratio favoring hypothesis SI over S2- 

Now since- P (SI | D) + P(S2|D) must equal one, the final 
values for posterior probabilities are determined by their 
ratio. It is no longer necessary to determine the prob- 
ability of the datum (P(D)), nor is it necessary to assess 
the probabilities P(D|S1) and P(DJS2) if the likelihood 
ratio I(D|S1:S2) is used. However, it may be easier to 
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compute the likelihood ratio based on the probability 
P (D | S) , which will be demonstrated in an example. 

The following example using this Bayesian method is 
presented. Three hypotheses under consideration are: 

51 = attack in Sector 1 

5 2 = attack in Sector 2 

53 = attack in -Sector 3 

The first difficulty encountered using this approach is 
determining prior probabilit ies. Should they be calculated 
using the first piece of evidence, or be based on other 
knowledge. In this example there is no reason to favor one 
over the other, so assign equal probabilities: P(S1) = P { S 2) 
= P (S3) = 0.333. 

The first piece of evidence received (D1) indicates 
strong possibility of attack in Sector 3 (S3) . To ease 
computation, all likelihood ratios will now be based on 
Sector 3, i. e. : L(D1|S1:S3), L(D1|S2:S5) and L(D1|S3:S3) 

will be used. Table 3 column L (D1 1 S (i) | : S3) lists the like- 
lihood ratios based on this first piece of evidence. 

The likelihood ratios can be computed in the following 
manner. Datum 1 (D1) indicated a high probability of attack 

in Sector 3, say p(D1)S3) = 0.8. To compute likelihood 

ratios, a comparison must be made with the other hypotheses. 
If D1 indicates a small probability of attack in Sector 1, 
(SI), say 0.10, then the likelihood ratio L(D1|S1:S3) = 

.1/.8 = .125. As previously stated, this likelihood ratio 
could have also been determined by saying that the prob- 
ability of attack in Sector 3 based on the datum is eight 
times greater than attack in Sector 1. 

The ratio for Sector 2 and 3 will be .2/. 8 = .250, or 

attack in Sector 3 is four times more likely than attack in 
Sector 2, based. on the datum. 
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The next piece of evidence received (D2) assigns the 
following likelihood ratios: 

L(D2|S1:S3) = .70/. 50 = 1 .4 
L (D2 | S2 : S3) = .80/. 40 = 1 .6 
L (D2 | S3: S3 ) = 1.0 

Now update all posterior probabilities for comparison 
using the three equations and three unknowns: 

P (S (i) j D) /P(S3| D) = [P (S(i) )/P (S3) ] • L(D1 |S(i) ) : S3) 

• L(D2 1 S (i) :S3) , for all i 

and solve simultaneously. 

Column P(S(i))10) of Table 3 (where D represents all 
evidence received) shows the final posterior probabilities. 
Now comparisons of hypotheses may be made using tnese poste- 
rior pr cbab ilities . 



TABLE 3 

BAYESIAH INFERENCE EXAMPLE USING LIKELIHOOD RATIOS 





P(S(i)) 


L ( D 1 j S (i) : S3 ) 


L (D2 | S (i) : S3) 


P(S(i) jD) 


S 1 


. 333 


. 125 


1.4 


.111 


S 2 


.333 


.250 


1 .6 


. 254 


S3 


. 333 


1.0 


1 . 0 


. 635 



P (S ( i) ) = Prior probability 

L (D 1 j 5 (i ) : S3 ) = Likelihood Ratio from rirst datum 

L (D2 j S (i ): S3 ) = Likelihood Ratio from second datum | 

P(S(i)|D) = Posterior probability | 

i 

i 
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There are several problems encountered using this 
approach. First, how should the prior probabilities be 
chosen? If they are equally likely, then equal probabili- 
ties could be assigned as done in this example. It may be 
better to use tne first piece of evidence receive, but these 
priors must sum to one. If the initial evidence indicates a 
strong possibility of attack in Sector 1, but lends no 
support to the other hypotheses, how should prior probabili- 
ties be assigned to the other hypotheses? Some prior prob- 
ability must be assigned to the other hypotheses, if not, 
they will always have a posterior probability of zero. 

Second, there can be a lack of consistency in assignment 
of likelihoods. There is no constraint to the magnitude of 
the ratio. In the example, a iikelihood ratio of 1.4 was 
used. Yet, if a strong indicator of attack (.9) was 
compared to a hypothesis with a very slight chance of attack 
(.001), then a likelihood ratio of . 9/. 001 = 900 would 

occur. This lack of constraint on the magnitude of the 
likelihood ratio may lead to inconsistencies as the infer- 
ence progresses through many likelihood ratios. 

Shafer and Tverskv £Bef. 13] remark that traditional 
Bayesian theory has been concerned with what they call 
observation design. This design deals with outcomes of 
statistical experiments. In the experimental space, the 
analyst knows the possible outcomes and answers. Prior 
probanilit ies for parameters can be assessed in advance. 
Bayesians have gradually extended their experimental space 
to the area of data analysis where probabilities are not so 
clearly defined. This possible over extension of Bayesian 
Inference could lead to its partial demise in the battle- 
field hypothesis space. 

This section has not been presented as a critique of the 
Bayesian method, but as an insight that there are problems 
with Bayesian methods as with Demps ter-Shaf er . 
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IV. A COHPUTATIOSAL VIES OF DEMPST ER-SHAFEE 

A major drawback to the use of Dempster-Shaf er has been 
the long calculation time required due to its computational 
complexity. For example, the computation of Belief requires 
time exponential in jTHETAj [Ref. 12]. As part of this 
thesis research a decision aid was created to assist in the 
Dempste r-Shaf er computations. Appendix A contains a 

description of the aid, as well as, a listing of the PASCAL 
code. This aid should not be considered a fully operational 
military decision aid, but rather a prototype or example of 
an automated evidence combination technique. 

The aid is designed to lead the user through the steps 
necessary to set up and use the De mpster-Shaf er theory. The 
aid's output has been verified by comparison with manually- 
computed solutions of problems using Dempster-Sh af er . 
Identical results were achieved. The output of the model 
used for these checks was Belief, Plausibility and Hass . 

Although this exponential computational time factor of 
Dempster- Shafer has been discussed in length, 

[Refs. 6,12,15], no actual computa tional data was found to 
support it. An additional benefit of the aid was the 
ability to new record these computational times. Also 
discovered through the use of the aid, was the memory limi- 
tation of the computer after all subsets of the Frame of 
Discernment were enumerated. 

Reduction in the computational complexity of 
Dempster-Shaf er will be addressed in this chapter, but 
first, a brief discussion of the De mps ter-Shafer aid created 
will be presented. Efficiency of the aid in terms of compu- 
tational time and memory requirements will be described in 
more depth. Then methods of reducing computational 
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complexity presented by Barnett [Ref. 12], and Gordon and 
Shortliffe, [Ref. 15], will be discussed. A third method of 
reducing computational time. Multiple Frames of Discernment, 
will be presented. 

A. A DEMPSTEB-SHAFER DECISION AID 
1. Des cr iption Of The Aid 

The aid is designed for a user familiar with the 
Dempster-Shaf er theory. The user will first see a screen 
informing him of the theory used in the aid. Next, the size 
of the Frame of Discernment is requested. A letter of the 
alphabet is assigned to each member of the Frame of 

Discernment. If a Frame had five members, THETA would be 
represented by (A,B,C,D,E) , or 7 members THETA = 

( A, B,C, E, E, F, G) etc. The user is then told the Frame will 
be represented by these letters. 

All subsets of the Frame are generated by procedure 
"Generate". For each subset, this procedure creates a 
PASCAL record structure that contains the items of Belief, 
Mass, ana Plausibility may be stored. The sets are sorted 
by size by procedure "Quicksort" to assist in the search 
efficiency throughout the program. Next; the aid informs 
the user of the item number for data input. The item number 
is just a means or keeping track of the numoer or loops 
through the program, which equals the number of data items 
combined. 

The user then enters the set for wnich he wants to 
assign Mass. Masses are assigned for all desired sets and 
then combined by procedure "Combine" using the 
Dempster-Shaf er theory. 

The user is then asked if he desires Beliefs to be 
computed. If so, 3elief and Plausibility are computed by 

procedure "Belief" and displayed. The program then returns 
to the input mode and will cycle until ended by the user. 
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Ef f iciency 0 £ The Aid 

As previously stated, a drawback to the use of 
Dempster-Shaf er is the time required to compute Belief. As 
discussed in Chapter III, Belief for a set requires the 
summing of the Masses of all subsets of the parent set. 
Figure 4.1 shows execution time for the computation of 
Belief by the aid, as well as, times for the other major 
program procedures. 

The execution time for procedure Generate increases 
exponentially. The procedure recursively generates subsets 
and must check for repetition of sets as it proceeds. The 
number of subsets is exponential in THETA, 2 exp |THETA|. 
As the size of THETA grows, the number of subsets grows 
exponentially, and also, the number of repeat subsets to be 
checked and eliminated increases. Fortunately, Generate is 
only executed once at the beginning of the session. It is 
possible to eliminate the generation of subsets altogether 
by reading in subsets from hard disk or floppy. All subsets 
for various sizes of Frames of Discernments could be stored 
and simply read at the beginning of the program. 

Unfortunately, the execution time pronlem for Belief 
computations is not so easily solved. Belief will be 
computed whenever the user desires a status of all the 
subsets. as the numner of subsets grows, this computation 
may cause the user an unacceptable waiting time to view 
Belief. During this time no input can De made for new 

Masses. Tne way this program is designed, multiple data 

entries can be made without involving the time intensive 
Belief procedure. In a battlefield scenario, the Belief 
computation would be done only when conclusions were 
required, not after each data input.- 

Perhaps a greater limitation in using micro 
computers to solve many real world problems is the large 
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PROCEDURE GENERATE PROCEDURE QUICKSORT 







Figure 4.1 



Execution Times Of Tne Aid. 
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NUMBER OF ITEMS IN FRAME OF DISCERNMENT n = NUMBER OF ITEMS IN FRAME OF DISCERNMENT 



memory requirements demanded by Dempster-Shaf er. With this 
particular coding scheme using a PASCAL record to represent 
each subset, memory limitations (320K RAM) were reached at a 
Frame of Discernment size of 9, which equals 2 exp 9 = 512 
subset of THETA. 

The military examples used so far in this thesis 
have had a Frame size of three elements and total size of 2 
exp 3=8. As seen in Figure 4.1, no noticable time delay 
occurs at this level, which is also well within memory 
storage limitations. However, the selection of Frame sizes 
above this level might cause unacceptable delays, especially 
in real-time battlefield applications. 

The three methods for reduction in complexity 
mentioned in the introduction will now be discused. 

B. BARNETT'S METHOD* SINGLETON HYPOTHESES AND THE IB 

COMPLEMENTS 

Barnett, [Ref. 12], showed that if all the subsets of 
the Frame of * Discernment can be reduced to singleton .hy poth- 
eses and their negations, computational time will be reduced 
from exponential to polynomial order. 

Before proceeding, a new military example will be used 
to demonstrate the formulation of Barnett's method. , The old 
example of sector of enemy attack is not weix suited for 
this metnod since it is difficult for the analyst to reduce 
the scope of the problem to singleton hypotheses and their 
comple ments. 

A suitable military example for the Barnett method is 
the analysis cf friendly axes of attack or enemy avenues of 
approach. Assume that the enemy or friendly force must 
choose the best approach into the combat area, given four 
avenues A, B, C and D. These singleton hypotheses and their 
complements will be the only sets considered. This example 
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will prohibit the analyst from assigning Mass to any combi- 
nation of elements, such as (A,D) , though he may assign Mass 
individually to (A) and (D) . Mass may be assigned in a 
sense to (B,C,D) by assigning Mass tc (not A) . 

If the problem at hand can be reduced to subsets using 
only tne hypotheses and their complements, this method does 
not restrict the analyst. He will see* reports about the 
avenues of approach into his sector or responsibility. 
These reports will consist of information about the terrain, 
weather conditions, natural and man-made obstacles, etc. 
The analyst will then assign Mass to each hypotheses based 
on the information and infer the most likely enemy approach 
or best friendly axis. 

1. Formulation Of Barnett* s Method 

Three steps will be used to represent the formula- 
tion of Barnett’s method: 

a) Combination of elements and complements. 

i) Combine evidence confirming each singleton 

hypothesis. 

ii) Combine evidence disconf ir ming each singleton 
hypothesis. 

Step 1 results in the formation of 2*N Belief 
functions (N= (THETA! ) , one for each element and 
one for each element complement. 

b) Combination of element pairs. Combine the confirming 
and disconf irming evidence for each element. This 
step forms N Eelief functions, one for each element of 
THETA. 

c) Combination within THETA. Combine ail elements of 
THETA tc produce one Belief function. 

Each step of this method will now be discussed. 
Only three avenues of approach (A,B,C) will now be consid- 
ered to simplify the computations. 
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a. Step 1: Combination Of Elements And Complements 

For each singleton, and for each complement, the 
Mass assigned by all evidence received is combined. Using 
avenue of approach (A) from the example, the evidence 
confirming (A) after receipt of two pieces of evidence would 
be: 



M * (A ) = Ml (A) a) M2 (A) = 1 - ( 1-MI ( A) ) • ( 1-H2 (A) ) 

(where (•) represents Mass at the end of Step 1) 

Since only evidence confirming (A) is used here. Mass is 
assigned by Ml or M2 only to (A) or to (THETA), see Figure 
4.2. No Mass will be assigned to (Null Set). ( 1 — H 1(A) ) and 
( 1- M2 (A ) ) represent Mass going to THETA from each individual 
assignment. Their product, (1-MI (A) ) •( 1-M2(A) ) represents 
combined Mass going to THETA. The remainder of Mass, 
1- (1-MI (A) ) • ( 1-M2 ( A) ) , goes to (A). No normalization is 
required since nc Mass goes to the Null Set. 

Suppose two more pieces of evidence are 
received, M3 and M4 , that disconfirm A. Mass from step 1 
assigned to (not A) would be: 

M * (not A) = 1- ( 1-M3 (not A ) ) • ( 1 -M 4 (not A)) 

1 - M* (not A) = Mass assigned to THETA 

This process continues for each element in 
THETA. The result is 2»N Belief functions, where each 
Belief function has two components one for each element and 
for each complement of the element. 

h. Step 2: Combination of Element Pairs 

In Step 2, the Mass for each element and its 
complement are combined. Figure 4. 3- shows an orthogonal sura 
example of the combination for Mass that supports an item 
and its complement. 
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Figure 4.2 Orthogonal Combination For Step 1. 
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Figure 4„ 3 Orthogonal Combination For Step 2- 
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Id this step Mass will be assigned to the 
element, complement of the element, THETA and the Nall Set. 
The Mass going to the Null Set: 

M n (N ull Set) = M* (A) «M * (not A) 

must be removed and so the Masses are normalized by the 
factor: 

1— M * (A) «M (not A) (1 - Mass going to the Null Set) 

Mass going to (A) in this step, M"(A), will be the combined 
Mass of the intersection of M* (A) and 1-M' (not A) , the Mass 
of the complement that was assigned to THETA from Step 1. 
The same analogy is true for the combined Mass going to (not 
A) in Step 2, M" (not A), see Figure 4.3. 

Ncw the number of Belief functions has been 
reduced from 2*N to N, one for each element and complement 
pair. The resultant Masses are represented by Barnett as p 
and c, where; 

p = M" (A) -= M* (A) ♦ ( 1-M ' (not A) ) / (1-M' (A) *M* (not A) ) 
c = M" (not A) = M* (not A) • ( 1-M* (A) ) / (1-M* (A) *M ' (not A) ) 

Also, the Mass assigned to THETA after combination in Step 2 
is represented by r: 

r = 1 - p - c 

and the variable, d, used later in Step 3 for normalization 
is the Mass assigned to the complement of the element and to 
THETA. "d" is used in Step 3 to compute the Plausibility of 
each element and is identified now for later computations. 

d = r + c 
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c. Step 3: Combination Within THETA 

Now the N Belief functions must be reduced to 
one by combination within the Frame of Discernment. The 
example Frame in this chapter: (A,B,C), will be used to 

demonstrate computations in Step 3. All results from Step 2 
will be represented as a function of their element, i.e.: 

p (A) = M" (A) 

c (A) = M M ( no t A) 

r (A) = 1-p <A)-c (A) 

d (A) = r (A) + c (A) 

The normalizing factor, K, for the overall 
Belief function is: 

K = 1 / (PRO D d (i) ) « ( 1 + SUM p (i) /d <i) ) (4.1) 

- EEOD c (l) , i in THETA 



or: 



K = 1/[d (A) *d (B) *d (C) ♦ (1 + p(A)/d(A) + p(B)/d(3) 

+ p (C) /d (C) ) - c(A)*c(B)*c(C) ] 

Belief for each item (i) is then computed by: 

Bel (i) = K ( p (i) • PROD d ( j) (4.2) 

+ r(i)* PROD c(j)), j#i 



then : 

Bel ( A) = K«p (A) *d(B)*d (C) + r (A) *c (B) *c (C) 

Belief for the complement of the item is computed by: 

Bel (not (i ) ) = K*[ (PROD d(i)*(SUM p(j)/d(j)) (4.3) 

+ c(i) *PRODd (j) - PROD c ( i) ], j*i 
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Bel (not A) = K[ (d (A) • d (B) «d (C) ) • (p ( B) /d (B) + p(C)/d(C)) 

+ c(A) •d(E) *d(C) - c (A ) *c (B) • c (C ) ] 

All Belief is computed in this manner to produce 
one Belief function. Evidential Intervals are then computed 
as follows: 

Evidential Interval (A) = [ Bel (A) , Bel (not A) ] 

The following numerical solution of the entire 
method using the example will be presented to aid in the 
understanding of this method - 

Eight items of evidence are received and Masses 
are assigned as follows: 
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0. 285 



r (A) =1 - p (A) - c (A) 

r (A) = 1 - 0. 429 - 0.286 = 

d (A) = r (A) + c (A) 

d (A) = 0.285 + 0.286 = 0. 571 

p(B) = (0. 8) • (1-0. 3) / ( 1- (0.8) • (0.3) ) =0.737 
c (B) = (0. 3) • (1-0. 8) / ( 1- (0.8) • (0.3) ) = 0.079 
r (B) = 1 - 0.737 - 0. 079 = 0. 184 
d (B) = 0. 184 + 0.079 = 0. 263 

p(C) = (0.838)* (1-0.2) /(I - (0. 838J • (0. 2) ) = 0.805 
c (C) = (0.2) • (1-0.838) /(I- (0. 838) • (0. 2) ) =0.039 
r (C) = 1 - 0.805 - 0. 039 = 0. 156 
d (C) = 0. 156 + 0.039 = 0. 195 

Step 3: 

K = 1/[ (. 571) * (.263) • (. 1 95) • (1 + . 429/. 571 + .737/. 263 

+ . 805/. 195) - (.286) •(. 079) *(. 039) ]= 8.998 

Computation of Beliefs: 

Bel ( A) = 3 .988* (.429*. 263*. 195 + . 285*. 079*. 039) = . 091 
Eel ( E) = 3 .988* (.737*. 571 *. 195 + . 1 84*. 286 •; 156) = . 360 
Bel(C) = 3.988* (.805*. 571 *.263 + . 1 56* . 2 86 * . 07 9) = . 496 
Bel ( not A) = 3. 988 (.029* (. 737/. 263 + .805/. 195) 

+ .286*. 263*. 195 - . 001 ) = .856 
Bel ( not B) =3.988* (.029* (. 429/. 571 + .805/. 195) 

♦ .014) = .620 

Bel ( not C) = 3. 988* (. 029* (.429/. 571 + .737/. 263) 

+ .014) = .467 



The Evidential Intervals would then be: 
A: T 0.09 1, 1-0.856] =[ 0.091,0.144 ] 

B: [ 0.360 , 0.380 ] 

C: [ 0. 496, 0. 533 ] 
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2. Computational Time For Barnett's Method 

A program was written to assist in the computations 
of Barnett's method. The program also enabled comparison of 
computational time and memory requirements with the normal 
Dempster-Shaf er method. A program listing for Barnett's 

method is contained in Appendix B. This program was written 
as a research tool and does not have the user friendly 
enhancements of the Dempster- Shafer decision aid in 
Appendix A. 

Figure 4.4 shows the timed computations for the 
Barnett method. The computer memory (320K RAM) could store 
a Frame of Discernment of up to 1000 items compared to 9 for 
full De mpst er-Shater . Figure 4.5 compares the Belief calcu- 
lation times for Barnett with the times already shown for 
the Dempster- Shafer aid. 

Barnett’s method offers a very appealing and time 
efficient use of Dempster- Shafer in a system where the 
following criteria are met: 

a) The Frame of Discernment can be adequately represented 
by the singleton elements and their complements. 

b) All evidence can be divided into confirming and 
d iscon f ir ming categories for each hypothesis. 

The example of the avenues of approach scenario 
could be well represented by this method. As shown, 
evidence based on terrain, weather, obstacles, etc., could 
be combined using Barnett's method. The resultant 
Evidential Intervals would then be compared to determine the 
most likely enemy route, or the best friendly route. 
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TIME IN SECONDS TIME IN SECONDS 






STEPS 1 AND 2: COMBINE WITHIN SET 




STEP 3: COMBINE WITHIN FRAME, 
COMPUTE BELIEF . • 




I 



Figure 4.4 Computation Times For Barnett's Method. 
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r 






COMPARISON OF BELIEF COMPUTATIONS 




Figure 4.5 Comparison Of Barnett Hith Full Dempster-Shaf er . 

C- GORDON AND SHORTLIFFE H ETHOD , HIERARCHICAL HYPOTHESIS 

SPACE 

Gordon and Shortliffe £ Ref. 15] agree that if evidence 
confirms or disconfirms singleton hypotheses, then Barnett's 
method produces the desired time reduction. However, there 
are many classes of problems where more flexibility is 
required. If the hypothesis space can be reduced to a 
strict hierarchy, many more real world problems could be 
handled. Gordon and Shortliffe use this Hierarchical 
approach in their work with MYCIN, a medical diagnostic aid. 
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but it is also well suited to a military example where 
intentions are being determined. 

The example of enemy attack can be expanded so tha 
analyst has tc consider the overall enemy intention: at 
reinforce, defend, delay, or withdraw. These intention 
be structured into a hierarchical tree such as tha 
Figure 4.6. Grouped in this manner, the example now 
the method proposed by Gordon and Shortliffe. 
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Figure . 4. 6 Enemy Intention Hierarchical Tree. 
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For mulation Of Gordon and S hor tlif f e ' s Method 

This hierarchical approach proceeds similarly to the 
singleton hypotheses method. All evidence is divided into 
confirming or disconf irming parts, only now there are pairs, 
triplets, etc. to which this evidence may be assigned. As 
in Barnett’s method, three steps will be used in the combi- 
nation cf evidence: 

a) Combine evidence for each element and each element 
complem ent . 

i) Combine all confirming evidence. Same as 
Barnett's first step. 

ii) Combine all disconfirming evidence 

b) Combine ail confirming evidence in THETA. 

c) Combine disconfirming evidence with confirming 
evidence from Step 2. 

Before proceeding, the hierarchical tree of Figure 
4.6 will be split into a tree with the elements of THETA, T, 
and one with the complements of T, Tc, see Figure 4.7 Since 
it is the superset of both trees, (THETA) itself, is not 
included in either tree. 

T = (AB, CDS, D£,A,B,C, D,E) 

Tc = (notAB, notCDE,notDE, notA,no tB, note, notD, notE) 

Evidence will be combined for |T| + JTcj items = 16, and 

final belief will be computed after Step 3 for |T| + THETA = 

17 it e ns versus 2 exp 5 = 3 2 items for the full 

Dempst er-Sha f er method. These reductions will decrease 

computational times and storage requirements. 

The steps for Gordon and Shortliffe's Method are as 

follows . 
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Tc 





Figure 4.7 Hierarchical Tree Of Elements And Complements. 

a. Step 1: Combine Evidence For Elements And For 

Complements 

Use Barnett' s equations to compute combined 
evidence for each element and each complement: 

M' (i) = 1- (1-MI (i) ) ( 1-M 2 (i ) ) , i in THETA (4.4) 

Where (') indicates Mass after Step 1 and (") will indicate 
Mass after Step 2. 

M ' ( A ) = 1 - ( 1 - Ml (A) ) • (1- M2 (A) ) 
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Each element (i) will now have a Mass M * (i) and a Mass 
1-M'(i) assigned to THETA . It is necessary to identify the 
Mass that each element’s combination sends to THETA for use 
in lat.er steps. Let M'i (THETA) equal the Mass assigned to 
THETA by element (i) during step 1. 

M» A (THETA) = 1-M' (A) 

The Masses for (not A) would then be: 

M* (not A) = 1-(1-M1(not A)) • (1 -M2 (not A)) 

M'notA (THETA) = 1-M* (not A) 

There are now 2«N Belief functions, N in ! and N 
in Tc. Step 2 will reduce the N Belief functions in I to 
one Belief function. 

b. Step 2: Combine All Confirming Evidence in 

THETA. 

All confirming evidence in THETA, all of which 
is in T , will now be combined. 

The combined evidence for THETA, M" (THETA) is: 



n” (THETA) = 



IT . 



r*r> r\ r\ 

c l\Kju 



M'i (THETA), i in T 



(4.5) 



This is the product of all Mass assigned to THETA by all 
elements in T. Now compute the combined evidence for all 
other elements in T: 



M" (i) = K*M»(i)*PROD M»j (THETA), i in T, (4.6) 

j not a superset of i 



The Mass of element (i) in Step 2 is- the product of its Mass 
with the Mass assigned to THETA by ail otiier sets in T 
except (i) , and supersets of (i) . Element (i) cannot 
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•combine with Mass that it assigned to THETA. In addition, 
it cannot receive Mass from a superset's assignment to 
THETA. The evidence that was assigned to a superset of (i) 
was not conclusive enough to assign it to element (i) . 
Allowing (i) to combine with M'j (THETA), (j) a superset of 
(i) , would allow assignment of Mass from which (i) was 
earlier restricted. This is not acceptable, under 
Dempster-Shaf er . The Mass for element (A) in the example 

would be: 

a" (A) = K* M' (A) •M' B (THETA ) «M' CD E ( THETA) • 

M'C(THETA) *M'DE (THETA) »M'D (THETA) •M'E (THETA) 

These calculations continue for all elements in 
T. The normalization factor, X, is the inverse of the sum 
of all the new Masses: 

K = 1/SUM M"(i), i in T or i = THETA (4.7) 

There are now N+ 1 Belief functions, N in Tc and 
one in T. Step 3 will reduce tnese to one Belief function. 

c. Step 3: Combine Disconf irming And Confirming 

Evidence 

Step 3 will now combine T and Tc to produce one 
Belief function. Each element of Tc will combine with T. 
Step 3 uses an approximation to 3, the orthogonal sum, which 
will be designated as a)'. M"(i) a)' M 11 (not i) will have non 

zero value on only (T union THETA), 

( ABCDE , AB, C D£ ,D£ , A , B , C , D, E) . Any Belief normally assigned 

by a) to g, g not in T, will instead, by 3', be assigned to 
the first ancestor of q in T. For example: 

(not E) a)* (CDE) = (ABCD) 3 (CDE ) 
but (ABCD) Z- (CDE) = (CD) , which is not in T. 
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In this case the Mass that would have been assigned to (CD) 
will be assigned to (CDE) the first ancestor in T. 

There are three cases for step 3 combinations. 
For each case there are two elements involved, one from Tc, 
(not Y) , and one from T, (X) . The combination case is based 
on the relationship of (Y) and (X). 

1. Case 1: X is a subset of Y: 

M (X ) = K«M" (X) •M' , notY (THETA) 

For example: 

M 11 (A) 3* M" (not AB) , A is a subset of AB 

M(A) = K«M" (A) *M"notAB (THETA) 

2. Case 2: X & Y = (Null Set) , i.e. X & notY = X 

a) Case 2(a): If (X union Y) in (T union THETA): 

M (X ) = K[M" (X) + M" (X union Y) •M' 1 (not Y) ] 

For example: 

K" d)' M" (not DE) , C & DE = (Null Set), 

CDE in (T union THETA), therefore: 

M(C) = K»[M’ ! (C) + M" (CDE) *M M (not D£) ] 

b) Case 2 (b) : If (X union Y) not in (T union THETA) 

then the Mass of (X) is not changed by the combi- 
na ticn : 

M (X) = K*M (X) 

M" (C) o)' M" (not E) , C & E = Null Set, 

CE not in (T union THETA), therefore: 

M (C) = K«M" (C) 

3. Case 3: X is a proper super set of Y 

a) Case 3(a): If X & nctY is a set in T: 

M(X) = K ♦ M" ( X) • M “not Y ( THETA) 
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M" (CDE) a)' M"(notC), CDE & ABDE = DE, a set in 
(T union THETA) 

M (CDE) = K«M" (CDE) •MnotC (THETA) 

b) Case 3 (b) : If X & notY is not a set in T (this 

case assigns Mass to the superset, which in all 
cases is X due to the strict hierarchy and unique 
parent requirements). 

M" (notE) a)' H" (CDE) , ABCD & CDE = CD, not in (T 
union THETA), but CDE is the first ancestor of CD 

M (CDE) = K* M (CDE) 

This process continues with each element from Tc combining 
with each element in T using one of the case rules. K, the 
normalizing factor is computed after each iteration of (not 
Y) from Tc. 

Normalization is done by summing the Masses in 
(T union THETA) and dividing all Masses by that value. 

K = 1/SUM M (i) , i in (T union THETA) (4.8) 



A shortcoming of the method described by Gordon 
and Shortliffe is that step 3 assigns all Mass to (T union 
THETA). No Mass remains in the- co mplement sets, therefore 
it is not possible to compute Evidential Intervals, 
[Bel(A) ,1-Bel(not A) ]. All comparisons between hypotheses 
must be done cn Belief alone. 

The following example of this method is provided 
for clarification. The hierarchical tree shown earlier in 
Figure 4.4 still applies. 

Masses are assigned as follows: 



M ( AB ) = 0.4, 

M (not CDE) = 0. 3, 
ii (DE ) = 0.6, 

M (A) =0.4, 



M AB (THETA) =0.6 
M note DE (THETA) =0.7 
a DE (THETA) = 0. 4 
M A (THETA) = 0.6 
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M(C) = 0.2, M C (THETA) = 0.8 

H(A) = 0.3, fl A (THETA) = 0.7 

There is only one element in Tc = (not CDE) , there are 
in T. 



Step 1: 

(A) is the only element with multiple Masses: 

M* (A) = 1- (ill (A) )* (1-M2(A) ) 

M*(A) = 1- (1-0. 4) •( 1-0.3) = 0.58 

M' A (THETA) = 1-0.58 = 0.42 

Step 2: Combine all confirming evidence in T. 

M" (T HETA ) = K *M * ABC (THETA) **M * D E ( THE TA) 

• M * A (THETA) «M *C (THETA) 

= K* (0. 9)* (0.5) • (0.42) • (0.8) = K*(0.081) 

M " (A E) = K «M ' (AB) *M * A ( THE TA) • M ' D E (THETA) • M ' C (THETA) 

= K« (0. 1 )• (0.42) • (0.5) • (0.8) = K*(0*.017) 

M M ( D E) = K>M ' (DE) *W A (THETA) • M'AB (THETA) 

= K» (0. 5) • (0. 42) • (0 .9) = K» (0.189) 

M M ( A ) = K« M* (A) • M * DE (THET A) 

= K* (0. 58) • (0.5) = K • (0 . 29) 

M"(C) = K*M* (C) • H* AB(THETA) »M * A(THETA) •M'DE (THETA) 

= K* (0.2)* (0. 9)* (0.42)* (0.5) = K*(.038) 

K = 1/The sum of all Masses in (T union THETA) 

= 1/a« (THETA) + M " ( A B) + M" (DE) + H"(A) + M"(C) 

= 1/(0.151 ) + (0.017) + (0. 189) + (0. 29) + (0. 038) 
= 1/0. 685 = 1.46 

Normalize the masses in T. 



four 
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The elements in T are: 



M" 


(THETA 


) = 


(1. 


46) • (. 151) 


= .221 


M" 


(AB) = 


(1 


. 46) 


• (.017) = 


. 025 


M" 


(DE) = 


(1 


. 46) 


• (.189) = 


. 276 


M” 


(A) = 


<1. 


46) • 


(.29) = 


.423 


M" 


(C) = 


(1. 


46)* 


(.038) = 


.055 



Total =1.0 



The elements in Tc are: 

M" (not CDE) = 0.3 M" note DE (THETA) = 0.7 



Step 

1 



3 : 



M" (notCDE) 


3' M" (THETA), 


THETA a super 


set 


of 


AB => 


Case 3 












THETA S 


(not CDE) = T 


HETA & AB = A3 


in 


(T 


union 


THETA) 












M (THETA) 


= K # M"(THETA) • M"not CDE (THETA) 


in 


this 


case, Mass 


would normally 


go to (THETA) 


and 


AB, 


see 


Fi-gure 4.2, 


here (THETA) 


its combined Mass 


for 


this 


iteration. 


and later when M" (not CDE) 




a 


« (AE) 



occurs, AB will receive its combined Mass. This is 
an iterative process versus the normal one step 
Dempster-Shaf er combination. 



a (THETA) = (. 221)* (.7) = 0. 155 (normalization will 

occur later) 



2. M" (not CDE) a)' H« (AB) , ABSCDE = Null Set, (AB union 
CDE) in (T union THETA) => Case 2 

M (AB) = [ M (AB) + M (ABCDE) •M (not CDE)] 

= [ .025 + (. 221) • (.3) ] = .091 

3. H" (not CDE) 3* H" (DE) , DE sunset of CDE => Case 1 

M (DE) = M" (DE) • MnotCDE (THETA) 

= ( . 276) • ( .7) = .193 
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4 



CDE&A 



M (not CDE) a)' M"(A), CDE&A = Null Set => Case 2, 
union CDE) not in (T union THETA) 

M ( A ) = M (A) (but will be normalized) 

= .423 

5. M" (not CDE) a)' M" (C) , C a subset of CDE => Case 1 
M (C) = M" (C) •H" notCDE (THE TA) 

= ( . 055) • ( . 7) = .039 

Compute K, the normalization factor. 

K = 1/ sum of Hasses in (T union THETA) 

= .155 + . 091 ♦ . 193 + . 423 + . 03 9 = . 901 



(A 



Normalized Hasses: 

M (THETA) = . 155/. 901 = . 1 72 
H (AB) = . 091 /.901 = . 1 0 1 

H (DE) = . 1 93/. 901 = . 2 14 

M (A ) = . 423/. 901 = .4 70 

M (C) = . 039/. 901 = .043 

Total = 1.0 

This process of combining elements from Tc with those in I 
would continue until all elements have been combined, in 
this example there was only one. Comparison of hypotheses 
would then be done based on final Belief values: 



Eel ( AB) 
Bel (DE) 
Bel (A) = 
Bel ( C) = 
Bel (B) = 



= M (AB) + M (A) = 

= M ( D E) = .214 
. 47 0 
. 04 3 

Bel (D) = Bel ( E) 



101 + . 470 =. 57 1 



= 0.0 



This would show strongest Belief (.571) in an aggressive 
enemy action (AB), which stands for attack or reinforce. 
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2. Computational Time For The Gordon and S ho rt li ffe 

Method 

A program was written to assist in the computations 
using the Gordon and Shortliffe technigue. A listing of the 
program is contained in Appendix B. This program, like the 
program for Barnett's method, does not have the user 
friendly enhancements of the Dempster-Shaf er decision aid in 
Appendix A. The purpose of the program was to assess the 
time calculation advantage of this method and determine 
memory reguir ements. 

Figure 4.8 shows a comparison of the Belief computa- 
tions for the three methods discussed so lar. Belief compu- 
tations for Gordon and Shortliffe are the same as the full 
Dempster-Shaf er, except there are fewer sets for which 
computations are necessary. • The time reduction occurs since 
many of the subsets of THETA are not considered. 

The maximum number of. subsets that can occur under 
the unigue parent restriction of the method are: (2«h - 1). 

This number results from a descending creation path from 
THETA, separating one element at a time. Now instead of 
having 2 exp N elements for which Belief must be computed, 
there are only (2»N-1). For a THETA of 10 elements, instead 
of having 2 exp 10 = 1024 elements, there may only be a 

maximum of 2*10 — 1=19 elements. This size reduction 

helps explain the computation times in Figure 4.8. 

The computer memory (320K RAM) could store up to a 
Frame of Discernment of up to 500 items using tne Gordon and 
Shortliffe method versus 1000 for Barnett and 9 for the full 
Dempster-Shaf er . 

Gordon and Shortliffe' s method is a practical and 
efficient use of Dempster-Shaf er in a system where the 
following criteria are met: 

a) A strict hierarchy of elements exists and each element 
has only cne parent. 
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COMPARISON OF BELIEF COMPUTATIONS 




Figure '4. 8 Comparison Of Belief Computational Times. 

b) All evidence can be divided into confirming and 
d isconf irming categories for each hypothesis. 

The example of the overall enemy intention fits well 
within the limits of this method. However, the comparison 

of the hypotheses of the intentions would have to be made on 
the result of Belief alone. The analyst would not know the 

Ignorance, 1-Eel (not X), remaining about each of the 

hypothesis. 

Shortcomings of this method a're: 
a) Loss of ability to compute Evidential Intervals as 
discussed earlier. 
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b) order dependence on combination. Gordon and 

Shortiifxe mention that their approximation to the 
Dernps ter-Shafer combination can be order dependent 
when a set and its parent have only one aescendent. 



D. MULTIPLE FRAMES OF DISCERNMENT 

One final method of reducing complexity, multiple Frames 
of Discernment, will now be discussed. This method may use 
any of the computational techniques discussed so far. Its 
reduction in complexity comes from the separation of the 
Frame into smaller, more manageable categories with which 
the analyst can work. Computation of Belief for these 
smaller Frames should now fall into the reasonable area of 
Dempster-Shaf er calculation times. 

The obvious requirement for this method is a logical 
separation of elements in the original Frame of Discernment. 
Once the elements are separated into multiple Frames, items 
from different Frames cannot be compared for they are now 
part of different Belief functions. 

As an example or this method, consider an expanded 
version of the enemy intentions problem. The analyst still 
desires to determine the overall intention. If the inten- 
tion is Attack, he wants to know which sector is most 
likely. He also desires to know if the enemy intends to use 
nuclear weapons, no matter what the tactic. The Frame of 
Discernment new has 9 items: 

1. Attack, Sector 1 

2. Attack, Sector 2 

3. Attack, Sector 3 

4. Reinforce 

5. Defend 

6. Delay 

7. Withdraw 
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but as described in 



8. Use Nuclear Weapons 

9. Not Use Nuclear Weapon 

A hierarchy of elements exists, 

Cnapter III with the sector example, unique parents do not 
exist. For example, let the Frame of Discernment be repre- 
sented as: A ,B ,C,D,E r F,G,H, I where: 

A = Attack, Sector 1 
3 = Attack, Sector 2 
C = Attack, Sector 3 
D = Reinforce 
E = Defend 
F = Delay 
G = Withdraw 
H = Use Nuclear Weapons 
I = Not Use Nuclear Weapons 

Now (ABC) is a multiple unique parent for (AB) , (BC) and 
(AC), rut (A) has two parents (AC) and (AE) . The same situ- 
ation exists for (B) and (C) . Therefore, Gordon and 
Shorliffe's method would not work here. 

The full Dempster-Shaf er method would create 2 exp 9 
subsets, which according to Figure 4.1 would require about 
100 seconds for Belief computations. Some of the pairs of 
elements of THETA, as discussed earlier in the Gordon and 
Shortliffe method, would not be of interest to the analyst. 
A solution to this dilemma is separating the problem into 
three Frames of Discernment. Each Frame will use one of the 
computational methods discussed. Any of the methods may be 
used where applicable in the multiple frames. All three are 
used here to show the diversity of this method. 

The overall intention of the enemy is still desired, so 
use the Gordon and Shortliffe method for the Frame: 

( AT , D, E ,F, G) , where (AT) is the attack intention and D 

through G remain as described above. 
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Use Barnett’s method for the question of the enemy’s 
intent to use nuclear weapons, (H,I), and use the full 
Dempster-Shaf er method to determine sector of attack. 
Caution must be exercised when discussing the most likely 
sector of attack for that Belief is conditioned on the 
Belief that the overall intent is attack. 

The timed calculations would now be: 



Full Dempster-Shaf er: 3 elements in THETA = approximately 

1 second 

Gordon and Shortliffe: 5 items in THETA = approximately 1 
second 

Barnett: 2 itmes in THETA = approximately 1 second 
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mplex and inefficient problem has now been reduced to 

manageable calculation for the intelligence analyst 
emps t e r- 5 haf er . 

deuce received while using this method of multiple 
cf Discernment, does not need to be separated into 
the multiple Frames. Cn the contrary, a report indi- 
an attack with nuclear weapons could be used for all 
mes if applicable. The Belief values of the various 
will net be compared and are not calculated using the 
of the other Frames. 

ermediate results can be saved on a disk and recalled 
re relevant data arrives. Therefore, one machine 
eep ail three methods running in an almost simulta- 
tate . 

method of Multiple Frames of Discernment is a viable 
tive to the full Demp ster-Shaf er method if the Frame 
rable into distinct categories. 
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V. SITUATION DEVELOPMENT AN ALYSIS. AN A PPL ICATION AREA FOR 

DEMPSTEB-SHAFER 



In Chapter II, the battlefield intelligence process was 
modelled as the Irt-erpr eta tion node in the command and 
control cycle. In this node, information was processed into 
intelligence by the analyst. Chapter III proposed the 
Dempster-Shaf er theory of evidence combination as an aid for 
the analyst in the Interpretation node. A decision aid and 
three computational techniques for reducing calculation time 
for the Dempster-Shaf er theory were presented in Chapter IV. 

Enhancement of a specific job application in the 
Intelligence node, through the use of the De mpster-S haf er 
theory and a decision aid similar to the one m Chapter IV, 
will now be analyzed. This specific job area. Army Division 
Situation Development Analyst (DSDA), was chosen due to its 
relevance to research conducted by MAJ L. Baltezore, U.S. 
Array, in his thesis at the Naval Postgraduate School 
[Ref. 16 ]. 

Baltezore's thesis proposed a Decision Support System 
hardware layout to assist the analyst in conducting situ- 
ation assessment at division level. A Knowledge Based 
System (KBS) was designed to conduct automated analysis 
concerning possible courses cf enemy attack. The use of 
Dempste r-Shaf er in this KBS will now be explored. Nhile not 
evaluating any specific technique, Baltezore proposed that 
some method of inference should be used to aggregate battle- 
field information stored in the data base. The Knowledge 
Eased System of the Division Situation Development Analyst 
and an intention assessment capability using the speed of 
Dempster-Shaf er . with small Frames of Discernment (Quick 
Assessment Capability), will be used as specific examples of 
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-actual Interpretation node duties that can be enhanced by 
Dempster-Shaf er . 

The jobs cf the DSDA will be briefly discussed as will 
the structure of the KBS. The use of a De mpster-S haf er 
theory in the Knowledge Based System as well as a Quick 
Assessment Capability through a decision aid will be 
evaluated. 

A. DIVISION SITUATION DEVELOPMENT ANALYSIS 

Situation Analysis at the division level is performed 
within the All Source Production Section (ASPS) of the All 
Source Analysis System (ASAS) , see Figure 5.1, [Ref. 16: 
p.25]« The ASAS is an Army project to establish a system 
that maximizes the productivity of Intelligence and 
Electronic Warfare (IEw). Through enhanced productivity, 
the Interpretation node will process information into intel- 
ligence more guickly, speeding up the command and control 
cycle. 

The Division Situation Development Analyst is usually 
the senior intelligence analyst in the All Source Production 
Section. Using a broad view of the enemy forces, the DSDA 
must determine key enemy objectives, rank potential enemy 
courses of action, and identify key targets, command and 
control nodes, or events indicating a specific course of 
action. The analyst bases his assessment of the enemy situ- 
ation on data passed through his work station [Ref. 16: p. 

32]. 

As described in Chapter III, tnis data is inherently 
uncertain and incomplete. Using this data, the analyst must 
interpret as much as possible about the enemy intent as 
quickly as possible. Figure 5.2- depicts the analyst's 
production cycle. The analyst would use the proposed 
Knowledge Based System to assist him in determining the 
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possible courses of enemy action. Within this structure 
exists the potential use of Dempster-S hafer . 
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Figure 5.2 Analyst's Production Cycle. 



B. THE KNOWLEDGE BASED SYSTEH 

As previously stated, the Knowledge Based System is 
designed to provide automated analysis concerning possible 
courses of enemy action. Analysis of the KBS will be 
limited to the structure of the system and the internal 
theory which correlates the data. 
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1. Knowledge Based S yst em St ru cture 

Baltezore’s proposed structure for the KBS is shown 
in Figure 5.3. This KBS interacts with the analyst through 
the Interrogation Module to determine the type of analysis 
desired. For example, is the enemy intention Attack, (one 
of the examples used in Chapter IV)? The Situation 
Assessment Processor accesses the Knowledge Base to deter- 
mine a rule list associated with the type of analysis 
desired. Using all data available, the KBS then establishes 
a probability value associated with the specified course of 
action. The Explanation Generator presents the analyst with 
the course of action considered, its probability value, and 
a rule audit trail of the deduction process. 



2. Knowle dg e B as ed Sy st em Theory 



The KBS deduction theory is based on work done by 
Ben-Basset and Freely, [Ref. 17], who proposed the use of 
classes, features, and relevancy pointers to conduct situ- 
ation assessment. 



a. Classes 

Classes are used to define battlefield situ- 
ations of interest to the analyst. For example, attack, 
reinforce, defend, delay, withdraw are enemy situations that 
may be represented as classes. If the general class is 
known, such as attack, then the specific location, such as 
the sector of attack, is desired through analysis. This 
characteristic is analogous to the Dempster-Shaf er frame of 
discern ment . 



b. Features 

Features are bits’ of information, such as the 
context of a report, related to the situation and used to 
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Knowledge Based System Structure. 
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An independent tank battalion is a 
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second echelon unit 
in a breakthrough 
classes is shewn in 



normally used as a lead attack element 
attempt. An example of features and 
Table 5. 4. 
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Class of Enemy Intention: Attack 
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Figure 5.4 Feature Probability Support For Class. 
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The significance of features and their support 
for the different classes under consideration are aggregated 
by some theory of evidence combination. Features, as infor- 
mation, may support several classes. The support of theses 
features for specific classes would now be combined to 
provide an overall indicator of support for the class. 

C. PROBLEMS WITH THE KNOWLEDGE BASED SYSTEM APPROACH 

The use cf a rule base and a priori probabilities for 
feature support of classes can be considered an expert 
system approach to the situation assessment task. The KBS 
conducts the actual analysis, the analyst merely gueries the 
system. The problem with an expert-type system approach is 
similar to that of a Bayesian approach, an endless task of 
defining and updating the rule base. 

Also, the rules used to determine the support of a 
feature for a class are not well defined. Table 5.4 showed 
hypothetical values of feature support for the attack inten- 
tion. The adjectives used, "massing" and "extensive", are 
very indistinct. The rules using them would then reguire 
some way of inferencing this indistinct adjective rrom the 
data case ("fuzzy sets"). For example, extensive prepara- 
tion would need to be defined as number of artillery rounds 
in an hour, cr number of targets engaged in a specific time 
period. If the data base supported these criteria, then the 
rules base could deduce that extensive preparation has 
occurred. 

Unfortunately, this use of explicit rules and indistinct 
features in an expert system can lead to a false sense of 
security of the battlefield. The analyst is dependent on 
the internal design of the system.. He is receiving the 
system's analysis ana may tend to discredit his own inter- 
pretation of the situation. Further, it is doubtful whether 
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any commander is presently ready to allow his analyst to 
accept a situation assessment from a "black box". The 
commander has too much risk associated with his decision 
based on this analysis to accept a fully automated analysis. 
Cushman , [Ref. 1], suggests that a commander will not be 
inclined to accept black box analysis , but instead support 
gradual automation based on transparent aids which allow the 
analyst to retain control. 

A decision aid, similar to that m Chapter IV, used in a 
quick assessment capability, is a complimentary alternative 
to the Knowledge Based System. This aid would act more as a 
parallel process with the analyst reflecting his view of the 
battlefield rather than an "expert" view given by the rules 
created by the "experts". The analyst is the expert in his 
division and cannot rely on rules or features created by 
other. 

Furthermore, the features suggested by Ben-Basset, 
[Ref. 17: p.486 ], were not intended for dogmatic application 
in all battlefield situations, but are given as a guide to 
the analyst for the analyst to use [Ref. 18: p. T-1]. The 
Dempster-Shaf er aid would be dependent on the analyst and 
his inputs rather than the analyst depending on the expert 
system. While the analyst would use these indicators as 
guides to assign Mass values, he would be free to change 
support values based on his knowledge of the situation and 
his prior experience. 

D. THE DEMPSTEB-SHAFER THEORY IN SITUATION DEVELOPMENT 

ANALYSIS 

The Dempster-Shaf er theory could be used in two ways in 
the DSD A. First, it could be used within the Knowledge 

Based Structure proposed by Baltezore. De mpster-Shaf er 
would be the combination technigue used to aggregate the 
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feature probabilities and produce the overall Belief for 
each class. The feature pr'obabilit ies would have to be 
represented as Masses and adhere to the definition of Mass 
presented in Chapter III. The use of the Gordon and 
Shortliffe method was discussed in Chapter IV using the 
example of determine overall enemy intentions. This tech- 
nique would reduce internal Knowledge Based System inference 
time and provide an efficient combination of evidence. 

Second, the aid could be used as a quick assessment 
device for the analyst. This use would be most relevant 
when the analyst can narrow the scope of the class. If the 
general class was attack and the analyst was concerned with 
the sector of attack, the D empster -Shafer aid would be used 
to determine Eelief for sector possibilities. This scope of 
use was represented through examples in Chapters III and IV 
using the sector of attack example. 

The aid would serve as a reflection of the analyst’s 
assessment of the battlefield as time progresses and reports 
are received. Using Belief values, the analyst will recog- 
nize the most likely sector of attack and advise the 
commander. Using Ignorance values, the analyst will reposi- 
tion or reorient sensors to investigate lack of knowledge of 
activit y. 

1. Adv an tages Of The Q.uick Assessment Use 

There are three main advantages to the quick assess- 
ment use over tne knowledge rased system: 

a) Absence of Rule Base 

b) Absence of Data Base 

c) Speed 



a. Absence of Rule .Base 

The analyst is not dependent on a predefined 
rule base to deduce support for the sector of attack. He 
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will consider each report and assign Mass values based on 
his- understanding of 'the enemy. This understanding will 
come from battlefield experience and knowledge gained from 
experts through instruction. 

b. Absence Of Data Base 

The reports received in the guick assessment 

capability will still be stored in some data base for future 
use. However, the analyst would not search through previous 
data, but interpret each piece as it is received. 

c. Speed 

For each .rule used, the KBS must satisfy the 
precedent {if portion) , to allow use of the antecedent (then 
portion). This process requires the continual search of the 
data base for conditions that satisfy the rule precedent. 
With a large rule case, such as that needed in the KBS 
structure to define all types of enemy activity, the cycling 
time would be prohibitive for a real time assessment 

capability. 

If time is available, the KBS procedure would be 
used to determine enemy intent, but if a quick assessment is 
desired, as discussed in Chapter II, the alternative capa- 
bility of Dempster-Shafer is the better option. 

E. RELEVANCY OF DEHPSTEB-S HAFER TO SITOATION DEVELOPMENT 
ANALYSIS 

It should be obvious that the examples used in preceding 
chapters are the same as the mission oojectives of the 
Situation .Development Analyst. Tne use of the 

Dempster-Shafer theory in this specific part cf the 

Intelligence node will speed the analytical process, 
providing the commander fast and accurate intelligence 
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This enhancement will reduce command and control 
and gain the benefit of this reduction discussed 
II. 



VI. SUHMABI AMD CONCLU SIONS 



A. SUHMABI 

This thesis has demonstrated the use of the 
Dempster-Sha f er theory of evidence in a decision aid to 
reduce command and control cycle time on the battlefield. 
The reduced cycle time allows the battlefield commander to 
interdict the enemy force earlier and gain a decisive 
advanta ge. 

The command and control cycle was modelled as a network 
to investigate data flow and network processes. The 
Interpretation node was determined to be the critical node 
and also a candidate for enhancement through application of 
the Dempster-Shafer theory. The Dempster-Shafer approach 
was presented as a plausible evidence combination technique 
when uncertain, incomplete, .and incorrect evidence must be 
comnined in a battlefield environment. 

Three methods for reducing the computational complexity, 
of the Dempster-Shafer theory, Barnett, Gordon and 
Shortliffe, and multiple Frames of Discernment, were demon- 
strated. These methods all have restrictions involving 
trade-offs between flexibility or scope and time efficiency. 
Military examples that met these restrictions were presented 
to demonstrate their possible use on the nattlefield. A 
decision aid based on the Dempster-Shafer theory was created 
and discussed. The aid eased the computational burden of 
Dempster-Shafer and allowed comparisons of computational 
speed with the three reduction methods. 

A specific application area for the De mpster-Shaf er 
theory. Situation Assessment in the All Source Analysis 
System (ASAS) • was described. The task of situation 
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assessment provided a good example of a process in the 
Interpretation node that can be enhanced by Dempster-Sh af er . 

B. CONCLUSIONS 

Dem pster- Shafer was not presented as a "cure-all" for 
evidence combination. Nor was it presented as a replacement 
for more pr cbabilisticly rigorous techniques such as, 
Bayesian Inference. The use of Dempst er-Shaf er in this 
thesis showed a logical method of combining data on the 
battlefield to help the analyst determine enemy intentions. 
The flexibility of Dempster- Shafer in handling battlefield- 
guality evidence should not be lost in discussions over its 
shortcomings. The intelligence analyst and battlefield 
commander need to make the best use of all data available in 
the most efficient and accurate manner possible. 
Dempster-Sh af er is a viable technique to assist in this 
process . 

The use of Dempster-Shaf er in the Knowledge Based System 
proposed by Baltezore would allow the intelligence analyst 
to conduct automated analysis. Tne analyst would access the 
data base of evidence and receive Belief values for his 
hypotheses of enemy intentions. The use of the Barnett, 
Gordon and Shortliffe, or multiple frames of Discernment 
methods, when applicable, would allcw for the most rapid 
computation of Belief. 

Automated analysis in this manner would allow access to 
mere data than the human processor could handle. Many 
analysts could use the same system over a period of time to 
analyze trends in the enemy activity. Furthermore, the 
analyst now has a backup system to his manual method of 
analysis. He can make the most use of the human-machine 
leverage discussed m Chapter II. 
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C. AREAS FOE FOBTHEB ANALYSIS 

This thesis has only began to investigate tne use of 
Dempster-Shaf er on the battlefield. The determination of 
mass values based on receipt cf evidence has been left to 
the analyst. The task of determining these values for 
Dempster-Shaf er, or any combination method, is formidable. 
Chapter V discussed the problem of a rule base with explic- 
itly assigned probabilities. Unfortunately there does not 
exist a data base from which probabilities of enemy inten- 
tions based on tactics can be extracted. 

Samet, [Bef. 19 ], has said that each sensor report has a 
reliability and accuracy associated with it. These features 
ccuid weight the ilass values assigned by the expert or 
analyst. The integration of these reliability and accuracy 
values into the Dempster-Shaf er mass values has not been 
discuss ed. 

The Knowledge Based System hardware of a Decision 
Support System, such as that proposed by Baltezore, must be 
designed to accept an evidence combination technique. 

These are but a few of the areas that are left for 
further exploration of De mps ter-Shaf er on the battlefield. 
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APP ENDIX A 

THE DE3PSTEB-SHAFEB DECISION AID 



This appendix contains the listing cf the code created 
for the decision aid and an explanation of the procedures. 
In general, the aid is a continuous loop that prompts the 
user for input from each knowledge source report received. 
It requests set identification and mass assigned to that 
particular set. The aid then combines the new input with 
the current mass values using Dempster-Shaf er . New Beliefs 
are computed only if directed by the user and then displayed 
to reveal the current status of knowledge about the Frame. 

The program is written in Turbo PASCAL for an IBM 

compatible computer. A main user routine and fourteen 
procedures, called throughout the routine,. make up the 

program. See figure A.1 for the program’s Flow Diagram. 

As stated, the program is composed of the following 

procedures. 

A. MAIN PROGEAM 

The main program, D5 ( Dem pster-Shaf er) , is an executive- 
like program that reacts to the user’s desires. Iz 

initially sets up the PASCAL Record Structure that will 
contain those items of information necessary for use in the 
program by Dempster-Shafer computations. These items are: 

1. Mass (Current Mass). 

2. lempmass (Temporary storage for Combined Mass). 

3. Newmass {User Mass input from Knowledge Source). 

4. Belief (Current Belief). 

5. ' ID (Set Describer) . 

6. Value (Number of Items in Set ID). 
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Figure A. 1 Flow Diagram for the Decision Aid Program. 
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An array is then created to allow for the maximum number 
of items in the Frame of Discernment which the computer can 
store. 

Next, the main program begins a calling sequence of the 
procedures. Within one of these procedures, many other 
procedures may be called tc accomplish the task at hand. 
After the aid is ready for user input, the main program asks 
if the user desires to continue. If the answer is yes, then 
the program cycles again. If no, then the user is allowed 
to request a final display of the output (Belief, 
Plausibility, and Mass) and the program terminates. 

The first procedure called by the program is Copyrite. 

B. PROCEDURE COPYRITE 

Procedure Copyrite is just a "cover sheet" for the aid 
identifying the theory used in the aid, and the creator. 

The next call is to procedure Initialize. 

C. PROCEDURE INITIALIZE 

Procedure Initialize reguests the number of items in the 
user’s Frame cf Discernment. It then determines the letter 
of the alphabet corresponding to the end of the size of the 
Frame, i.e., 5 = E, 7 = G. The letters of the alphabet will 
then be used as identifiers for each separate set. This use 
of letters allows for set operations within PASCAL. 

The number of sets that can be made from the Frame is 
then computed, where size = 2 exp N, N = number of items in 
the frame. 

Procedure Set Up is called next. 
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PROCEDURE SET UP 



This procedure sets up 
Frame of Discernment. Tho 
letters "A" through the end 
"Z" , used only fcr easy iden 
Procedure Levels is next 

E. PROCEDURE LEVELS 

Procedure Levels is used 
nations of set size N split 
1. In other words, K item 
time, through N items lata 
The results of tnis proc 
are used to help reduce "do 
where it is necessary to ch 
as when Belief is computed, 
determined by its Value (num 
the sets with value less th 
possibility. This action e 
parent sets that cannot be s 
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Figure A. 2 



Execution Times Versus Frame Size. 
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The computation time is high since identical sets are 
created in the recursive calls. Duplicate sets must be 
checked for and eliminated. As the frame size grows, this 
time factor is compounded, see Figure A. 2. 

The next procedure called is Quicksort. 

G. PROCEDURE QUICKSORT 

Quicksort, like Generate, is a recursive routine that 
sorts the sets according to the number of items in the set 
ID. This resulting sorted array of records allows for a 
more efficient check for subsets. A subset can only occur 
at a lower level of Value (number of items in set ID) than 
its superset. Therefore, loops computing Belief or 
Plausibility need only check lower levels of sets and not 
the entire tree of sets. 

Quicksort is almost unaffected by the number of items in 
the Frame of Discernment due to its efficiency, see Figure 

A. 2. 

If the user desires to continue with the program. 
Procedure Entermass is then called. 

H. PROCEDURE ENTERMASS 

Entermass is the main user input procedure that assigns 
new mass values to those sets identified by the current 
knowledge source. 

The first time called, Entermass allows the user to 
request a display of all the various sets that can occur. 
The user may then print this list for future reference. 

The aid tnen prompts the user for the identity of the 
set and the mass to be assigned to that set. Set identities 
are checked to insure they are in the domain of the Frame of 
Discernment. ilass is checked for containment in the set 
[ 0 , 1 ]. 
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Total mass entered during the current session is checked 
to make sure it does not exceed 1.0. Set ID are also 
checked to see if mass has already been assigned to the set 
during this session. If so, the user can change the mass 
value, or leave it at the current value. 

When the total mass entered eguals 1.0, the new masses 
are combined with the old mass value. 

This combination is done using Procedure Combine. 

I. PROCEDURE COMBINE 

This procedure combines masses using Dempster’s Rule of 
Combination and then normalizes the new masses. Current 
masses { mass es at the end of the previous session) are then 
replaced with the new combined mass values. 

Time of execution of the procedure is reduced by 
limiting the combination to only those sets that have a mass 
from the previous session, and those sets that were assigned 
masses during the current session. This efficient operation 
keeps execution time to a minimum, see Figure A. 2. 

The inter secti on-set of sets with masses are than iden- 
tified and assigned the new masses. After all orthogonal 
sums (See Chapter III), are computed, new masses are assigned 
back to tne Record for each set. 

If there is mass assigned to the null set, by default 
the mass assigned to the other sets does not total 1.0, then 
the masses are normalized. 

The Normalizing Factor is displayed to the user. New 
masses are calculated and assigned to the sets. The user is 
then asked if computation of Beliefs is desired. 

J. PROCEDURE BELIEF 

The computation of Belief is a very time consuming 
process since the mass for all subsets of each set must be 
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summed. The execution time for Procedure Belief is shown in 
Figure A. 2. 

Both Value (part of the Record) and Levels (created by a 
procedure) are used to reduce the search time for subsets to 
cnly those sets with the possibility of being a subset. 
Only those sets at lower levels of the Record structure may 
be subsets. After Belief values are computed, they are 
assigned back to the Record for each set. 

The Eelief, Plausibility, and Mass for each set are then 
displayed to the user. 

K. PROCEDURE DISPLAY 

Display is a procedure that writes ail sets, from higher 
order to lower, to the screen with their respective Eelief, 
Plausibility, and Mass. However, only those sets with 
Eelief values greater than zero are displayed. The screen 
will display only 12 sets at a time to allow easy viewing 
and, if desired, printing by the user. 

After this step, the masses from the next knowledge 
source are then entered and the process continues. 

L. ADDITION AI PROCEDURES 

There are several procedures used throughout the program 
not mentioned above that are described here. 

1. Procedure Display 2 

Display2 is the procedure used in the first itera- 
tion of Enter mass that, if requested by the user, displays 
all the sets (cnly set IDs). 

2. Procedure Che ckanswe r 

Checkanswer is a procedure used to insure the answer 
to Yes or N o guestions is in the correct form. 
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3. Procedure Max integer 

Maxinteger is used in the Combine Procedure to limit 
the search for the set containing the result of the inter- 
section of two masses. Since the sets are sorted, it is 
cnly necessary to look in the lowest level (or highest value 
in a high to low sorted array) of the two set IDs combined. 

4. Procedure Fa c tor ial 

This procedure is used in Procedure Levels ana 
computes the factorial of a number. This result is used in 
the computation of combinations for the values in the Levels 
array. 
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Listing of the PASCAL program used in the Decision Aid. 



program ds; 
const 

tempsize = 512; 



type 

node = record {stores data relative to d-s] 

mass : real; 

tempmass : real; 

newmass : real; 

belief : real; 

plausible : real; 

id : set of char; 

value; integer; 
end ; { node } 

branch = array[ 1 .. tempsize] of node; 



var 

count, ol.dcount , iteinnumber , numi terns , ii, size,cntr, one 
integer ; 

k, enditem : char; 
tree. : branch; 

goodanswer, continue' ; boolean; 

ch44 , contanswer , beliefanswer : char; 

l. evelend ; array[1..20] of integer; 

tempmasscntr : array[1..25] of integer; 

oldmasscntr : arrayCl .. tempsize] of integer; 

{include the following procedures m this main program} 

{ $ 1 be lief . pas } 

{$i factor ial . pas } 

{ $ i .l.eve l s . pas } 

{$i copyr i te . pas } 

{$i generate . pas } 

{$i setup. pas} 

{$i checkanswer . pas } 

{$i initialize . pas } 

{$i quicksort . pas } 

{$i sort. pas} 

{$i display . pas } 

{$i disp!ay2 . pas } ’ 

{$i maxmteger . pas } 

{$i combine. pas} 



lid 



{$i entermass . pas } 



begin {main} 

itemnumber : =0 ; {keeps track of the number of evidence 
combined} 

{begin calling sequence} 

copyrite ; 
initia lize ; 
setup; 

levels ( numitems ) ; 
ch44 ; = ' a ' ; 



write In; 
textcolor ( 14 ) ; 

write.!n( ' Generating Subsets '); 

generate ( tree[ 2 ] , ch44 ) ; 
wntel.n ; 

writel.n( ' Sorting Subsets '); 

textco lor (15); 



one : = 1 ; 

quicksort (one, size ) ; 

for li : = 1 to size do {initialize set values} 
begin 

tree[ii] . belief : = 0.0; 
tree [ ii ]. mass : = 0.0; 
end; { for } 

c Irscr ; 

wntel.n ; 

27 i tz . 1 n. ( 1 *************************************************** 
*★★***★★**• 

wri teln ( ' Your Frame of Discernment has been expanded into 
all Subsets ' ) ; 

write In ('and the Null Set'); 

wnte!n( 'As Each Item Of Evidence Is Received, You Will Be 
Prompted ' ) ; 

wnteln( 'For The Mass Distribution'); 

********** • ^ . 
writel.n; 

{check for users desire to continue} 
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continue: = true; 
while continue do 
begin 

repeat {until, goodanswer} 
write.l.n; 

write.ln (' ====> Do You Wish To Continue? Y or N'); 
textcol.or ( 12 ) ; 

writel.n(' A No Answer Will Exit Program'); 

textcol.or (15); 

read ( kbd , contanswer ) ; 

write.l.n (contanswer); 

checkanswer ( contanswer , goodanswer );{ checks to see 
if in Y,y,N,n) 

until goodanswer; 

if contanswer = 'n' then 
continue: = false; 

if continue then 

entermass ( itemnumber , numitems ) ; 
write.ln; 

repeat {until, goodanswer) 
write.l.n; 
wr itel.n ; 

wntel.n ('====> Do You Wish To Display Belief? Y or 

N ' ); 

read ( kbd , be l ief answer ) ; 

writel.n ( belief answer ) ; 

checkanswer ( be l lef answer , goodanswer ) ; 

until goodanswer; 

if be lief answer = ' y* then 
begin 
belief; 
display ; 
end ; { if } 

end; {while ) 

end . {main } 

procedure setup; {creates initial. 2 sets, null, and 
theta(entir frame)} 

begin 

tree[l].id := ['z']; {null set) 



112 



treeCl ] . value : = 0; 

tree[2].id := [ ' a ' . . enditem] ; {entire frame} 

tree[2 ]. value : = numitems; {number of elements in set id} 

cntr:= 2; 

end; { setup} 

procedure initialize; {user procedure to determine size of 
frame } 

var 

l, j : integer ; 
ok : boolean; 

begin 

for j:= 1 to 5 do 
write.l.n; 

repeat {until ok} 

writeln (' ====> Enter Number of Items in Frame of 
Discernment ' ) ; 

{$1-} readln (numitems) {$1+}; 
ok:= (I0resu.lt = 0); 
if not OK then 
begin 

text.col.or ( 12 ) ; 

writel.n( '*** Improper answer, retry'); 
textcolor ( 15 ) ; 
end ; { if } 

until ok; 

writeln (numitems); 

enditern:= chr ( numitems + ord('a')-l); 
writeln; 

write In ('***. Frame of Discernment has numitems :3, ' items 

and will be listed'); 
writeln ; 

wntel.n ('as Character Set a through ' , enditera ); 
size:= 2; 

for i:= 1 to numitems-1 do {computes 2**n} 
size:= size * 2; 

end; [initialize] 

procedure copynte; {cover display} 
var 

11 : integer;, 

dummy : char; • 
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begin 

textco.l.or ( 15 ) ; 
textbackground ( 1 ) ; 
cl.rscr ; 

for 1.1. :=1 to 5 do 
write.l.n; 
write.ln( ' 

CREATED TO COMBINE 
write.l.n; 
writel.n( ' 

THEORY OF EVIDENCE 
write.l.n ; 
write!n( ' 

' ) ; 



YOU ARE ABOUT TO USE A DECISION AID 



1 ); 



EVIDENCE USING THE DEMPSTER-SHAFER 



1 ); 



COMBINATION 



write l.n; 
write.l.n ; 
textco l.or (12); 

wnteln( ' CAUTION l l : Only those willing to venture 
beyond ' ) ; 

wntel.n( ' Bayesian Inference should 

continue ' ) ; 

textcolor ( 14 ) ; 
write.l.n ; 
write l.n? 

write!n( ' Written by CPT William H. Cl.eckner, US Army 
as a prototype ' ) ; 

wnteln( ' decision aid for combining intelligence and 
determining an.' ) ? 

write.l.n( ' enemy commanders main attack sector.'); 
for U.;= 1 to 5 do 
write In; 
textcolor (15); 

write.l.n (' Push any key to continue'); 

read (kbd, dummy); 

cl.rscr; 



end ; 

procedure generate(var rl : node; var beginitem : char); 
{conducts a recursive call to generate subsets} 



var 

i , j : integer ; 

chl,chll : char; 

newset : boolean; 

begin 

for chl:= beginitem to enditem do {remove a character 
at a time until null.set reached} 
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begin 

if chi in rl . id then 
begin 

tree[cntr + l].id:= rl . id - [chi]; 
if tree[cntr + l].id <> [] then 
begin 

cntr:= cntr + 1; 
newset:= true; 

for j : = 1 to cntr - 1 do {check for repeat 

character sets} 

begin 

if (tree[j].id = tree[cntr].id) and newset 

then 

begin 

newset:= false; 
cntr:= cntr - 1; 
end; { if } 
end; { for} 

if (newset) then {if sec not created before, 
generate its subsets} 

beg in 

tree[cntr] . value : = rl. value - 1; 
chll:= succ(chl); 
generate ( tree[cntr ] , chll); 
end ; {if} 

end; {if} 
end ; { if } 
end; { for } 
end; {generate} 

procedure quicksort ( var first ; integer; var last : 
integer ) ; 



var 

l, j, dividing.line : integer; 
temporary : node; 

begin 

i:= first; 
j:= last; 

dividing line ; = tree[ ( first + last) div 2]. value; 
repeat 

while tree[i], value > dividingline do 
i : = i + 1 ; 

while tree [ j ]. value < dividingline do 

•j:= j-i; 

if l <= j then 
begin 
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end; 



temporary := tree[i]; 
tree[i]:= tree[j]; 
treeCj]:= temporary; 
1 : = i+1; 
j:= j-1; 



until. 1 > j; 

if first < j then quicksort ( first , j); 

if i < last theh quicksort(i, last); 

end; {quicksort} 

procedure entermass ( var item : integer; var numitems : 

integer ) ; 

{enter new evidence temporary sets} 
label 

la, 20 ; 



var 

goodentry , goodanswer , ok ; boolean; 

disp.l.ayanswer , ch3 , ch33 , changeanswer : char; 

check, tmass : real; 

qq, nn, kk , j j , min, i : integer; 

tempset : set of char; 

tempid ; stnng[20]; 



begin 

item:= item + 1; 
c.l.rscr ; 

write.l.n ('*** You are entering the mass distribution for 
item number : ' ) ; 
write.l.n ; 

write In (' — > ' , item :4); 

write.l.n; 

if item = 1 then {ask user if needs display of 

sets } 

begin 

repeat {until goodanswer} 

wntel.n (’====> Do you need a display or all tne 
sets? Y or N ' ) ; 

write In; 

wntel.n (' This Feature Will Only Appear 

Once so Use ' ) ; 
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sets 



YyNn} 



writeln ( ' PrtSc key for hard copy of 

' ); 

writeln ; 

read ( kbd , displayanswer ) ; 
writeln (displayanswer); 

checkanswer (displayanswer , goodanswer ) ; {insure 



'Y') then 



until goodanswer; 
c.lrscr; 

if (displayanswer = 'y') or (displayanswer = 

begin 

display2; {display if requested} 
end; {if} 



end; {if} 



for nn:= 1 to size-1 do 
begin 

tree[nn] . tempmass : = 0.0; 
tree[nn] . newmass : = 0.0; 
end ; { for } 



repeat 



goodentry:= true; 
check := 0.0; 
count ;= 1; 

while check < 0.9999999 do {end while is label 20} 
begin 

tinass:= 0.0; 

repeat {until goodanswer label 10} 
tempset:= []; 
writeln; 

repeat {until ok} 

writeln ('====> Enter Set ID : '); 

{$1-} readln (ternpid) {$1+}; 
ok:= (iOresul.t = 0); 
if not OK then 
begin 

textcolor ( 12 ) ; 

writel.n( ' *** Improper answer, retry'); 
textcolor (15); 
end; { if } 
until, ok; 

writeln (tempid); 



117 



for ii:=l to length ( tempid ) do {check for 

set in frame) 

begin 

if (tempid[ii] in [ ' a ' . . enditem] ) and 
not (tempid[ii] = ' ') then 

goodanswer:= true 
else 

begin 

goodanswer:= false; 
textcol.or (12); 
write.l.n ; 

write.l.n ('*** set ID not in Frame, 

Retry ' ) ; 

textcol.or (15); 
goto 10; 
end ; { if } 

tempset:= tempset + [ tempid[ ii ] ] ; 
{create tempset of id) 

end ; { for } 

{check for attempt to enter mass for set already accessed) 

for qq:= 1 to size-1 do 
begin 

if (tempset = tree[qq].id) and 
( tree[qq] . tempmass>3 . 0e+00 ) then 

begin 

repeat {until goodanswer) 
write l.n ; 
textcol.or (12); 

writel.n(' You Have Already Entered A 
Mass For Set ' , tempid,' = ' , tree[qq] . tempmass :6:3); 

write.l.n ( ’ Do You Wish To Change This 

Input? Y or N ' ) ; 

textcol.or (15); 

read (kbd , chang eanswer ) ; 

checkanswer ( changeanswer , goodanswer ) ; 
until, goodanswer; 

if (changeanswer = ' y') or 
(changeanswer = 'Y') then 

begin 

check :=check - tree[qq ]. tempmass ; 
goto 10; 
end 
else 

goto 20; 

end; { if ) 
end ; { for } 
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10 : 



until goodanswer; 



repeat {until ok) 
write. l.n? 

write.ln ('====> Enter Kjss for Set'); 

{$1-} read.ln (tmass) {$1+}; 

ok:= (iOresult = 0); 
if not OK then 
begin 

textco.lor ( 12 ) ? 

writeln('*** Improper answer, retry'); 
textco.lor (15); 
end; { if } 

until ok; 

write.ln (tmass :8:5); 

for jj:= 1 to cntr-1 do {assign new masses to 

tree } 

begin 

if- tree[jj].id = ternpset then 
begin 

tree[ j j ] . tempmass : = tmass; 
tempmasscntr [count ]: = jj? 
end; { if } 
end; { for ) 
write.ln ; 

count : = count + 1; 

check := check + tmass; 
textcolor ( 14 ) ; 

write ('Total Mass = ', check :8:5); 
write (' (*** Reminder: Total Mass Must be = 1.0 

to exit loop ***)'); 

textcolor (15); 
write In; 

20 : end; {while) 

if (check <= U. 9999999) or (check >= 1.0000001) then 
begin 

goodentry:= false; 
write.ln; 
textcolor ( 12 ) ; 

write.ln ('*** Warning, Mass Total. = ', check :4:2,' 

Is Greater Than 1.0 ***'); 



119 



write.ln; 

writel.n ('*** You Must Input All Masses For This 
Item Of Evidence ***'); 

textco.lor ( 15 ) ; 
writein ; 

for i:= 1 to count do 

tree[ tempmasscntr [ i ] ] . tempmass: = 0.0; 

end; {if} 

until goodentry; 
if item > 1 then 

combine; {return newmasses in tree} 

if item = 1 then {do not combine} 

begin 

o!dcount:= count; 
for kk ; = 1 to count do 
begin 

tree [ tempmasscntr [kk ] ] . mass : = 
tree[ tempmasscntr[kk] ]. tempmass ; {reassign newmasses to 
tree } 

oldmasscntr [kk] : = tempmasscntr [kk ] ; 
end; { for } 
end; {if} 

end; {entermass} 
procedure combine; 

{combines evidence using dempster-shaf er } 
var 

h, hh , j , j j , i , li , k , kk , j2 , high , tempcount : integer; 

norma.l.f actor , totalmass, tempmass : real; 

intersection : set of char; 

tempoldmasscntr : arrayC 1 . . tempsize] of integer; 
begin 

tex tco lor ( 14 ) ; 

write l.n ( ' Combining Masses ' ) ; 

textcol.or (15); 

tempmass := 0.0; 
totalmass := 0.0; 
tempcount := 0; 

{ check for intersection and increment newmass by mass 
product } 



120 



for j:= 1 to oldcount do 
begin 

for i:= 1 to count do 
begin 

intersections tree[o.ldmdSscncr [ j ] ] . id * 
tree[ tempmasscntr [ i ] ] . id; 

• if intersection <> [] then 

begin 

max integer ( terapmasscntr[i] , o.l.dmasscntr [ j high ) ; 

for jj:= leve.l.end[tree[high] . va.lue+1 ]+l to 
leve.l.end[tree[high] . value] do 

begin 

if tree[jj].id = intersection then 
begin 

tree[ j j ] . newmas s ; = tree[ j j J . newmass + 
tree[o.ldmasscntr [ j ] ] . mass * tree[ tempmasscntr [ i ]]. tempmass ; 

tempcount:= tempcount + 1; 
tempoldmasscntr [ tempcount ] ; = j j ; 
for li ; = 1 to tempcount-1 do 

if tempoldmasscntrC ii] = jj then 
tempcount := tempcount-1; 

end; {if} 
end; { for } 
end; {if} 
end ; { for } 
end ; { for } 

oldcounts tempcount; 
for kk : = 1 to oldcount do 

ol.dmasscntrLkk] := tempo l.dmasscntr [kk J ; 

for j 2 : = 1 to size-1 do 
tree[ j2 ] . mass := 0.0; 

for h:= 1 to oldcount do {assign combined masses back to 
tree } 

begin 

totalmass:= totalmass + tree [o.l.dmasscntr [h] ]. newmass ; 
tree[ol.dmasscntr[h] ] .mass : = 
tree [o l.dmasscntr [ h] ] . newmass ; 
end ; { for } 

{if mass assigned to nullset (by default that totalmass 
assigned to } 

{sets <> l.O) then normalize] 

if totalmass < 1.0 then 
begin 

normal factor : = totalmass; 
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texrcol.or (14); 

write In (’ Normalizing Factor = normalf actor :G:3); 

textcol.or ( 15 ) ; 

for k:= 1 to o.l.dcount do 

tree[oldmasscntr[k] ] .mass : = 

• tree[oldmasscntr[k] ] . mass / norma.l.f actor ; 
end { if } 

else 

wntel.n( ' Normalizing Factor = ' ,totalmass :6:3); 
end ; { combine } 

procedure belief; {sums masses of all subsets and assigns to 
parent } 

var 

l , i2 , n , n2 , not va l.ue : integer; 

d : char; 

notset : set of char; 

begin 

textcol.or (14); 

wnte.l.n( ' Computing Beliefs '); 

textcol.or (15); 

for n:= 1 to size-1 do 
begin 

tree [ n] . be l.ief : = tree [n] . mass ; 

for n2:= (( levelend[ tree[n] . value] )+l ) to size-1 do 
{only checks its mass(n) and tnose sets with} 
begin 

if tree[n2].id <= tree[n].id then 

tree[ n] . belief : = tree[n] . belief + tree[ n2 ] . mass ; 
{if subset then increment belief} 
end ; { for } 
end ; { for } 

for l ; = 2 to size-1 do 
begin 

notset:= [ ' a ' . . enditem] - tree[i] . id; 
notva.lue;= abs ( tree[ i ] . va.l ue - numitems); 

for i2:= ( leve lend[ notva.l.ue + 1] + 1) to 
l.evelend[notvalue J do 
begin 

if tree[i2].id = notset then 

tree[ i ]. plausible : = 1 - tree[ i2 ] . bel ief ; 
end ; { for } 
end ; { for } 
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end ; { belief} 

procedure display; {displays set. ids, beliefs, masses} 
var 

jj : integer; 
ch, dummy : char; 

begin 

c.l.rscr ; 

for jj:= 1 to size-1 do 

if tree[ j j ]. belief > 0.0 then 
begin 

textcolor ( 13 ) ; 

wnte.ln ( ' 

■ ). 

textcolor ( 15 ) ; 
write ( ' Set ID; ' ) ; 
textcolor ( 14 ) ; 
for ch:= 'a' to enditem do 
begin 

if ch in tree[jj].id then 
write ( ch ) ; 
end; { for} 
textcolor ( 15 ) ; 
write ( ' Belief = 

' ) ; textcolor ( 14 ) ; write ( tree[ j j ] . belief : 6 : 3 ) ; 
textcolor ( 15 ) ; 
write (' Plausibility = 

' ); textcolor ( 14 ); write ( tree[ j j ]. plausible ; 6 ; 3 ) ; 
textcolor (15); 
write ( ' Mass = 

' ) ; textcolor ( 14 ) ; write ( tree[ j j ] . mass : 6 : 3 ) ; 
textcolor ( 15 ) ; 
write.I.n; 

if (jj mod 12) = 0 then {allows scrolling to stop 
temporarily } 

begin 

textbackground ( 14 ) ; 

write ('Push Any Key to Continue Display'); 
textbackground ( 1 ) ; 
write l.n; 

read ( kbd , dummy ) ; 
end; {if} 
end ; { for } 
end; {display } 

procedure disp.!ay2; {displays only set ids to choose from in 
procedure entermass} 

var 
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jj ; integer; 
ch ; char; 

begin 

c.l.rscr ; 

for jj:= 1 to size-1 do 
begin 

write ( ' Set ID; ' ) ; 
for ch:= 'a' to enditem do 
begin 

if ch m tree[jj].id then 
write ( ch ) ; 
end; { for } 

writel.n; 
end ; { for } 
end; { display 2 } 

procedure levels (var itm ; integer); 
var 

tleve.l. : array[1..20] of integer; 
i : integer; 

numerator : real; 

begin 

l.evel.endC itm + 1 ];= 0 ; 
numerators factorial ( itm) ; 
for i:= 1 to itm do 

tlevel[i];- trunc ( numerator / ( factorial ( i-1 ) * 

factorial ( itm-i+1 ) )); 

for i;= 1 to itm do 

levelend[ itm-i+1 ]: = tlevel[i] + levelend[ itm-i+2 ] ; 
end; { level end} 



procedure checkanswer ( var answer ; char; var test : 
boolean ) ; 

{checks to see if answer in set Y y N n returns false if 
not } 

begin 

if answer in ['y','Y','n', , N'] then 
tests true 
else 

beg in 

tests false; 

wnte.ln ( ' *** Answer was not in acceptable form, retry'); 
end; {else} 
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end; { checkanswer } 

procedure maxinteger ( var a : integer; var b : 

largest : integer); 

{determines the maximum of two integers (used 
combine ) } 

begin 

if a > b then 
largest := a 

else 

largest := b; 
end; {mininteger} 

function factorial (n : integer) : real; 
var 

fact,j : real; 
l : integer; 

begin 

fact:= 1; 

for i:= 1 to n do 

fact:= fact * i; 

factorial.:- fact; 

e nd ; 



integer; var 
in procedure 
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APPENDIX B 

THE METHODS OF BARNETT AND GORDON AND SHORTLIFFE 

This appendix contains the listing of the code created for the 
methods of Barnett and Gordon and Shortliffe described in Chapter 
IV. Note that Step 2 of Gordon and Shortliffe is not complete. 

It was written to compute Belief computation times. 
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Listing of the PASCAL program for Barnett’ 

program dsbarnett; 

const 

tempsize = 100; 



type 

node = record 

id : set of char; 

p : real; 

c : real.; 

r ; real.; 
d : real; 
ptemp : real; 

ctemp ; real.; 
bel : real; 

belcomp : real; 
end ; { node } 

branch = array [ 1 .. tempsize] of node; 

var 

position, item, numitems, size, i, j, k 
dprod, cprod, normfac, s, sprime, ratio 
real.; 

tree : brancn; 
tempid : stnng[20]; 

.finished, done : boolean; 

donans, f mans , n, y , enditem, dummy, ch : 

begin 

write.l.n (' enter number of items in frame’); 
read.l.n ( numitems ) ; 

enditem: = chr( numitems + ord('a')-l); 
size:= numitems; 

for i:= 1 to size do 
begin 

with tree[i] do 
begin 

id:= [cnr(i + ord ( ' a ' ) -1 ) ] ; 

p : = 0.0; 

C : = 0.0; 

r : = 0.0; 
d : = 0.0; 
ptemp := 0.0; 
ctemp:= 0.0; 



s method. 



: integer; 

, bigc, bigd 



char ; 
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end; {with} 



end; { for } 



done:= false; 
repeat { until done} 

finished:= false; 
repeat{until finished} 
writeln (' enter set id'); 
readln ( tempid ) ; 

for j: = 1 to size do 
begin 

if [tempid] = tree[j].id then 
positions j; 
end; { for } 

writeln (' enter mass that supports: tempid); 

readln ( tree [posit ion] . ptemp ) ; 

writeln (' enter mass that supports compliment of: 
readln ( tree[position] . ctemp ) ; 
writeln ; 

writeln (' Step One and Step Two '); 

with tree[position] do 
bey in 

s : = l-((l-ptemp) * (1-p)); 
sprime: =1- ( ( 1-ctemp) * (1-c)); 
p: = ( s* ( 1 -sprime ) ) / ( 1-s* sprime ) ; 
c : = ( sprime* ( 1-s ) ) / ( 1-s* sprime ) ; 
r : =l-p-c ; 
d: =c+r ; 

writeln ( 'p,c,r,d',p,c,r,d); 
end; {with } 



write In (' finished? y or n'); 
readln ( f mans ) ; 

if ( f inans = 'y') or (finans = 'Y') then 
finished:= true; 

until, finished; 
writeln ; 
writeln ; 

write In (' Step Three '); 



bigc: =1.0; 
bigd:= 1.0; 
ratio: =0.0; 



' , tempid ) 
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for k:= 1 to size do 
begin 

with tree[k] do 
begin 

ratio:= ratio + p/d; 
bigd:= bigd * d; 
bigc:= bigc * c; 
end; {with} 
end; { for } 

ratio:= 1 + ratio; 

normfac:= l/((bigd * ratio) -bigc ) ; 

write.l.n (' ratio, normfac, bigd, bigc, 

' , ratio, normfac, bigd, bigc) ; 
for i:= 1 to size do 
begin 

cprod ; = 1.0; 
dprod ; = 1.0; 
for j:=l to size do 
if j <> i then 
begin 

cprod := cprod * tree[j].c; 
dprod ;= dprod * tree[j].d ; 
end ; { if } 
with tree[i] do 
begin 

bel: = normfac* (( p*dprod ) + (r*cprod) ) ; 
belcomp:= normfac * ( bigd * (ratio - 1 

c*dprod - bigc); 

end; {with} 
end; { for } 
read In ( dummy ) ; 
c.lrscr; 

for j:= 1 to s.ze do 
begin 

textcolor (13); 

write l.n ( ' 

' ); 

textcolor ( 15 ) ; 
write (' Set ID: '); 

textcolor ( 14 ) ; 
for ch:= 'a* to enditem do 
begin 

if ch in tree[j].id then 
write ( ch ) ; 
end ; { for } 
textcolor ( 15 ) ; 
write ( ' Belief = 

' ) ; textcolor ( 14 ) ; write ( tree[ j ] . bel : 6 : 3 ) ; 
textcolor ( 15 ) ; 



- (p/d)) + 
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write (' Plausibility = ' ) ;textcolor(14) ;wnte(l- 
tree[ j ] . be.l.comp :6:3); 
wnte.ln; 

if (j mod 12) = 0 then {allows scrolling to stop 
temporarily} 

begin 

textbackground ( 14 ) ; 

write ('Push Any Key to Continue Display'); 

textbackground ( 1 ) ; 

wnte.ln; 

read (kbd, dummy ) ; 
end; { if } 
end ; { for } 

write l.n (' done-end program, y or n?'); 
read.ln ( donans ) ; 

if (donans = 'y') or (donans = 'Y') then 
done:= true; 

until done; 

end. {main } 



Listing of the PASCAL program for Gordon and Short.liffe 
method. 

program usshort; 
const 

tempsize = 100; 



type 

node = record 

id : set of cnar; 

p : real; 

c : real; 

rp : real; 

rc : real; 
ptemp : real; 

ctemp : real; 

be l : real; 

end ; { node } 

branch = array [ 1 .. tempsize] of node; 



var 
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‘ tempset ; set. of char; 

position, item, numitems, size, i, j, k ; integer; 
dprod, cprod, normfac, s, sprime, ratio, bigc, bigd 
real. ; 

tree : branch; 

setid, tempid : string[20]; 

finished, done ; boolean; 

donans, f inans , n, y , enditem, dummy, ch : char; 
begin 

write.l.n (' enter number of items in theta'); 
read.l.n ( numitems ) ; 

enditem:= chr (numitems + ord('a')-l); 

write.l.n (' enter number of sets to include theta'); 

read.l.n ( size ) ; 

for i:= 1 to size do 
begin 

with tree[i] do 
begin 
p : = d.0; 
c ; = 0.0; 
rp:= 0.0; 
rc ; = 0.0; 
ptemp:= 0.0; 
ctemp:= 0.0; 
be.l : = 0.0; 
end; {with} 
end ; { for } 

tree[l3.id;= [' a '.. enditem] ; 
for l ; = 2 to size do 
begin 

tempset :=[]; 

write.l.n (' enter set id for set number ',i); 
readl.n ( setid ) ; 

for j:= 1 to length ( setid ) do 

tempset ;= tempset + [setid[j]]; 

witn tree[i] do 
begin 

id:= tempset; 
end; {with } 

end; { for } 



done:= false; 
repeat {until done} 
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finished:= false; 
repeat. { until finished} 
write.l.n (' enter set id'); 
read l.n ( tempid ) ; 

tempset : =[ ] ; 

for j:= 1 to length ( tempid ) do 

tempset: = tempset + [tempidCj]]; 

for j:= 2 to size do 
begin 

if tempset = tree[j].id then 
positions j; 
end; { for } 

write l.n ( 1 enter mass that supports: ', tempid ) ; 

read l.n ( tree [posit ion] . ptemp ) ; 

write.l.n (' enter mass that supports compliment of: ',terupid) 

read.l.n( tree[position] . ctemp) ; 
write.l.n ; 

write l.n (' Step One '); 

with tree[position] do 
begin 

s : = 1- ( ( 1 -ptemp) * (1-p)); 

sprime : =1- ( ( 1-ctemp) * (1-c)); 

p: = ( s* ( 1 -sprime ) ) / ( 1-s* sprime ) ; 
c : =( sprime* ( 1-s ) ) / ( 1-s* sprime ) ; 
rp: =l-p ; 
rc : = 1 - c ; 
end ; [with } 



write l.n (' finished? y or n'); 
read l.n ( f mans ) ; • 

if (finans = 'y') or (fmans = 'Y') tnen 
finisned:= true; 

until finished; 
writel.n; 

write.l.n ( 'Step Two ' ) ; 

for i:=l to size do 
begin 

tree[l].p:= tree[l] ,p * tree[i].rp; 
write.l.n ( ' p, rp ' , tree[ 1 ] . p, tree[ i ] . rp ) ; 
end ; { for } 

for i:= 2 to size do 
beg in 

for j:= 2 to size do 
begin 
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if(j <> i) or not (tree[i] . id<=tree[ j] . id) then 
tree[i].p:= tree[i].p * tree[j].rp; 
end; { for } 
end ; { for } 

normfac:= 0.0; 

for i:= 1 to size do 

normfac:= normfac + tree[i].p; 

write-in ( ' normfac ' , normfac ) ; 
normfac := l/normfac; 

for i:= 1 to size do 

tree[i].p:= tree[i].p * normfac; 

writel.n ( 1 Step 3 '); 

for i: = 2 to size do 
begin 

for j:= 2 to size do 
begin 

if tree[j].id = ['a' ..enditemj - tree[i].id then 
tree[j].p:= tree[j].p * tree[i].c; 

if tree[j].id <= tree[i].id then 

tree[j].p:= tree[j].p * tree[i].rc; 
end ; { for } 
end ; { for } 

normfac : =0 . u ; 
for i:= 1 to size do 
begin 

normfac := normfac + tree[i].p; 
end ; { for } 

for j:= 1 to size do 
begin 

tree[ j ] . p : = tree[ j ] . p/ normfac ; 
end; { for } 

for i:= 1 to size do 
begin 

tree £ i J . be 1 : =0 . 0 ; 
for j : = 1 to size do 
begin 

if tree[j].id <= tree[ij.id tnen 

tree[ l] . be.l. : = tree[i].bel + tree[j].p; 
end ; { for } 
end ; { for } 
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c.l.rscr ; 

for j:= 1 to size do 
begin 

textcol.or ( 13 ) ; 

write In ( ' 

'); 

textcol.or (15); 
write ( ' Set ID: ' ) ; 

textcol.or ( 14 ) ; 
for ch:= 'a' to enditem do 
begin 

if ch in tree[j].id then 
write (ch ) ; 
end ; { for } 
textcol.or ( 15 ) ; 
write ( ' Belief = 

' ) ; textcol.or (14); write ( tree[ j ] . be l : 6 : 3 ) ; 
textcol.or (15); 

write (' Mass = '); textcol.or ( 14 ); write ( tree [ j ]. p 

: 6 : 3 ) ; 

write l.n; 

if (j mod 12) = 0 then {allows scrolling to stop 
temporarily } 

begin 

textbackground ( 14 ) ; 

write ('Push Any Key to Continue Display'); 

textbackground ( 1 ) ; 

write.l.n; 

read (kbd, dummy ) ; 
end; { if } 
end ; { for } 

write.l.n (' done-end program, y or n? ' ) ; 
read.l.n ( donans ) ; 

if (donans = 'y') or (donans = ' Y ' ) then 
done:= true; 

until done; 

end . {main } 
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